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ABSTRACT 


A  lobe  sweep  interferometer  system  is  described  which  has 
been  developed  at  the  Boeing  Scientific  Research  Laboratories. 
Continuous  measurements  are  made  of  the  intensity  of  the  total 
coherent  noise  received  from  the  sun  at  221  me  and  the  phase 
difference  between  the  signals  from  two  antennas.  A  phase  lock 
method  of  maintaining  a  1  kc  difference  between  two  mixer  injection 
signals  at  251  me  is  described.  A  phase  compensator  is  discussed 
in  which  a  fixed  phase  relationship  is  maintained  between  a 
reference  signal  and  the  1  kc  signal  resulting  from  a  fictitious 
source  at  the  predicted  position  of  the  sun.  It  is  shown  that 
the  use  of  the  compensator  simplifies  angle  measurements  of  burst 
positions  and  facilitates  measurement  of  the  position  of  the  radio 
center  of  the  sun.  Although  several  aspects  of  the  system  are 
not  unique,  most  of  the  system  is  discussed  for  the  sake  of  con¬ 
tinuity  and  completeness.  A  discussion  is  given  regarding  the 
interpretation  of  the  received  data  to  provide  position  measurements 
of  bursts  in  relation  to  solar  co-ordinates.  The  effects  of  fringe 
blurring  and  correction  for  an  extended  source  distribution  are 
discussed.  Several  factors  affecting  the  precision  of  position 
measurements  are  presented. 


PREFACE 


A  lobe  sweep  Interferometer  system  has  been  developed  at  the 
Boeing  Scientific  Research  Laboratory  to  be  used  primarily  for 
the  study  of  electromagnetic  emissions  originating  from  the  sun 
at  meter  wavelengths.  The  purpose  of  this  report  is  to  communicate 
to  the  reader  the  design  philosophy  used  in  the  development  of  the 
Boeing  system,  some  of  the  practical  limitations  of  the  system, 
and  how  the  recorded  data  are  interpreted.  It  is  intended  that 
this  report  will  form  the  framework  for  reports  to  follow  which 
will  deal  specifically  with  the  reception  and  interpretations 
of  solar  data. 

As  it  presently  operates,  the  lobe  sweep  interferometer 
records  continuously  the  intensities  of  _  ^lar  radio  bursts  at 
221  me  and  their  positions  relative  to  the  radio  center  of  gravity 
of  the  sun.  The  measurements  to  date,  which  have  been  made 
principally  to  establish  the  feasibility  of  the  system,  have 
indicated  the  desirability  of  further  refinements.  These  refine¬ 
ments,  which  are  scheduled  for  early  incorporation  into  the  system, 
include  increasing  the  baseline  length,  installing  preamplifiers 
directly  at  the  antenna  terminals,  and  the  design  and  construction 
of  a  device  which  modifies  the  operation  of  the  output  phase 
meter.  This  device,  called  a  phase  compensator,  is  currently 
undergoing  system  tests.  Because  of  its  uniqueness,  the  phase 
compensator  is  discussed  in  considerable  detail  in  the  report, 
even  though  it  has  not  yet  been  finally  incorporated  into  the 
lobe  sweep  system  for  routine  solar  observations. 
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PART  I 


HISTORICAL  SURVEY  OF  INTERFEROMETER  SYSTEMS 

Except  for  inherent  simplicity,  the  total  power  interferometer 
has  certain  disadvantages,  especially  at  low  signal- to-noise 
levels.  The  shortcomings  consist  primarily  of  poor  amplitude  and 
angular  resolution  in  the  presence  of  non-coherent  noise.  In  the 
measurement  of  the  position  of  short-term  radio  noise  enhancements 
of  the  sun,  Little  and  Payne-Scott  (l)  employed  a  lobe  sweep  inter¬ 
ferometer  in  which  the  lobe  pattern  from  a  two  element  interferometer 
was  swept  back  and  forth  across  a  source  at  a  constant  rate  by  the 
use  of  a  mechanical  phase  shifter  inserted  in  one  of  the  local 
oscillator  lines. 

An  improvement  in  the  capability  of  measuring  the  position  and 
intensity  of  weak  sources  was  made  possible  with  the  development 
of  a  new  type  of  interferometer  system  by  Ryle  (2).  Several 
configurations  are  possible  in  thiB  new  type  of  system,  but  all 
are  characterized  by  the  feature  that  the  output  recorder  mean 
deviation  is  independent  of  the  internal  noise  of  the  system.  The 
mean  deviation  depends  on  the  intensity  of  the  source  and  the 
path  difference  between  the  two  interferometer  antennas.  The 
statistical  fluctuations  about  the  mean  are  a  function  of  the 
total  noise  in  the  system.  These  systems  usually  have  two  basic 
components,  a  device  which  places  a  tag  on  the  coherent  noise 


signal  in  the  form  of  amplitude  modulation,  and  a  synchronous 
detector  which  compares  the  received  modulation  component  with 
a  reference  signal  of  the  same  frequency. 

The  phase  switch  system  developed  by  Ryle  (2)  U6ed  a  180  degree 
phase  switch  inserted  in  one  of  the  antenna  legs.  The  phase  switch 
varied  at  a  step  rate  and  produced  a  switch  frequency  modulation 
component  in  the  system.  The  modulation  component  was  synchronously 
detected  at  the  switch  frequency  producing  an  output  which  was  a 
measure  of  the  source  intensity  and  the  position  relative  to  the 
interferometer  lobe  pattern.  An  extension  of  the  phase  switch 
system  to  measure  phase  and  therefore  position  continuously  was 
made  by  Penfield  (3)  who  used  a  servo  phase  control  to  maintain 
a  null  at  the  output  of  the  synchronous  detector. 

Hanbury  Brown,  Palmer  and  Thompson  (4)  msed  the  lobe  sweep 
principle  in  addition  to  employing  a  synchronous  detector  to 
observe  the  position  and  intensity  of  weak  radio  stars.  Crystal 
controlled  local  oscillator  signals  were  used,  which  were  separated 
by  a  small  difference  frequency.  This  accomplished  the  same  end 
result  as  the  use  of  a  mechanical  phase  shifter  and  eliminated  the 
problems  associated  with  mechanical  limitations.  Cohen  (3)  used 
a  somewhat  similar  instrumentation  scheme  in  polarisation  measure¬ 
ments  at  200  me. 

A  further  improvement  in  the  lobe  sweep  system  can  be  realized 
by  using  local  oscillators  which  are  phase  locked  to  a  constant 
difference  frequency  as  was  done  with  the  Vanguard  Minitrack 
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system  (6,  7).  This  was  accomplished  in  the  Minitrack  system  by 
using  reactance  control  in  a  voltage  tunable  L/C  oscillator  circuit. 

Continuous  phase  measurements  using  a  lobe  sweep  system  have 
been  made  by  Fleischer  (8)  who  used  a  phase  meter  in  place  of  a 
synchronous  detector.  This  instrument  is  used  to  study  the  position 
and  motion  of  long  duration  bursts  at  51 7  me. 

A  swept  frequency  technique  discussed  by  Wild  and  Sheridan  (9) 
has  been  used  to  measure  the  position  and  sise  of  bursts  originating 
from  the  sun  over  the  range  from  40-70  me.  The  operating  principle 
relies  on  the  fact  that  the  path  difference  between  the  receiving 
antennas  is  a  function  of  the  wavelength  which  is  swept  over  a  band 
at  a  fixed  rate. 

Multi-element  interferometers  have  been  used  to  measure  directly 
the  brightness  distribution  across  the  solar  disk  at  meter  wave¬ 
lengths  by  Boischot  (10)  and  Firor  (11).  In  the  multi-element 
interferometer  system  the  sun  moves  across  narrow  strips  which 
remain  stationary  in  the  sky  and  result  from  the  antenna  polar 
pattern.  Since  the  angular  separation  between  these  strips  is 
greater  than  the  angular  size  of  the  sun,  the  system  can  un¬ 
ambiguously  resolve  several  centers  of  activity  emitting  simul¬ 
taneously,  provided  the  secondary  lobes  are  sufficiently  attenuated. 
The  main  shortcoming  of  the  multi-element  interferometer  is  in  the 
resolution  of  short  duration  bursts.  Since  the  scan  time  is  long 
compared  to  the  duration  of  short  term  bursts,  no  accurate  position 
measurements  of  single  transient  events  can  be  carried  out.  For 
example,  an  intense  burst  occurring  at  a  position  of  a  secondary 
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lobe  response  will  appear  to  an  observer  to  be  a  low  intensity 
burst  at  a  position  corresponding  to  the  position  of  the  main 
lobe.  The  position  of  noise  storms  can  be  determined  indirectly 
by  making  a  statistical  analysis  of  a  number  of  separate  scans. 

Two-element  interferometers  have  been  used  to  resolve  short 
duration  bursts  by  Owren  (12)  and  Fokker  Cl3 )  using  total  power 
and  phase  switch  interferometers  respectively.  Both  observers 
determined  the  position  of  active  centers  by  noting  the  times 
and  hence  the  position  when  the  enhanced  emission  component 
was  either  in  or  out  of  phase,  or  at  phase  quadrature.  Owren 
measured  the  positions  relative  to  the  position  of  the  quiet 
sun  which  was  assumed  to  coincide  with  the  optical  center.  In 
later  measurements  Fokker  measured  the  position  relative  to  the 
optical  center  of  the  sun  as  distinct  from  the  radio  center  of 
gravity. 


PART  II 


DISCUSSION  OF  THE  BOEING  SYSTEM 


General  description  of  the  system 

The  principle  components  of  the  Boeing  system  are  shown  In 
the  simplified  block  diagram  given  in  Figure  2.1.  The  lobe 
sweep  results  from  the  fact  that  the  two  local  oscillator  in¬ 
jection  signals  at  approximately  251  me  are  separated  by  1  kc. 

Since  phase  lock  is  used,  the  difference  frequency  stability  is 
determined  by  the  1  kc  reed  oscillator  shown  in  the  figure.  An 
AGC  system  is  employed  which  operates  from  the  coherent  noise 
signal  output  from  a  narrow  band  1  kc  filter.  Intensity  recordings 
are  made  over  the  range  from  the  level  of  the  quiet  sun  to  power 
increases  of  35  db  above  this  level.  This  range  is  covered  by 
two  recorders,  one  for  low  level  signals  and  one  of  extended 
range.  A  third  recorded  output,  resembling  in  form  a  sine  wave 
and  herein  called  the  interferometer  output,  is  that  from  the 
synchronous  detei  or.  This  output  is  used  chiefly  as  an  indication 
of  proper  system  operation. 

Hie  fourth  recorded  output  is  from  the  phase  meter  and  is 
an  analog  representation  of  the  phase  difference  between  the 
1  kc  coherent  noise  signal  and  the  1  kc  reference  signal  from 
the  phase-locked  oscillators.  The  recorder  trace  is  basically 
a  sawtooth  waveform.  The  phase  compensator  shown  in  the  figure 
is  currently  being  tested  and  has  not  yet  been  incorporated 
into  the  system  for  routine  solar  observations.  It  represents 
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an  improvement  over  the  present  phase  meter  system,  in  that  is 
will  remove  the  phase  discontinuity  from  the  present  sawtooth 
waveform  and  it  will  enable  the  phasemeter  to  indicate  the 
variation  in  position  of  the  radio  center  of  gravity  of  the  sun. 

The  phase  compensator  accomplishes  these  by  inserting  into  the 
system  a  time  varying  phase  which  is  nearly  equal  to  that  predicted 
for  the  center  of  the  sun  as  it  moves  relative  to  the  antenna 
baseline.  As  the  phase  compensator  is  considered  to  be  unique, 
it  is  discussed  in  detail  in  Part  IV. 

In  the  event  of  a  radio  burst  the  1  kc  signal  amplitude 
and  phase  will  be  the  resultant  of  the  superposition  of  the 
signal  from  the  quiet  sun  (pre-burst  level)  and  that  of  the 
burst.  By  measuring  the  total  intensity  and  phase  of  these 
superposed  signals  the  intensity  and  phase  of  the  burst  are 
calculated.  A  discussion  later  in  the  report  will  be  concerned 
with  the  interpretation  of  the  data  to  provide  burst  positions 
and  intensities. 

A  summary  of  the  opei ating  characteristics  and  the  parameters 
of  the  complete  system  is  given  in  Table  2.1. 

Site  location 

The  receiving  site  is  located  in  a  thinly  populated  area 
approximately  11  air  miles  from  the  Boeing  Scientific  Research 
Laboratories  in  a  location  comparatively  removed  from  power 
lines  and  roads.  Early  field  tests  had  established  that  the 
level  of  radio  noise  interference  at  meter  wavelengths  was 
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TABLE  2.1 


SUMMARY  OF  OPERATING  CHARACTERISTICS 
FOR  THE  BOEING  LOBE  SWEEP  SYSTEM 


Nominal  operating  frequency 
Lobe  sweep  frequency 
IF  frequency 
System  noise  figure 
Over-all  receiver  bandwidth 
Antennas 
Gain 

Separation 

Regularly  recorded  outputs 
Intensity 

Extended  range  Intensity 
Analog  phase 

Available  recording  time 
constants 


221.54  me 
1  kc 
30  me 

4.5  -  5.0  db 
640  kc  or  125  kc 

Two  10-element  Yagjs ^equatorial  mounts 

Nominally  12  db  above  an  isotropic 
radiator 

109.0  feet  ■  33 .2  meters 

80°/o  full  scale  »  100  x  10‘22wnf  2 (cps)"1 
(quasi-logarithmic ) 

Full  scale  >  31,600  x  l('"22wm*2(cps)"1 
(quasi-logarithmic ) 

Proportional  to  the  phase  difference 
between  the  two  receiving  antennas 

.25,  .5,  1.0,  2,  and  4  seconds 
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sufficiently  low  that  the  site  would  be  satisfactory  for  solar 
observations. 

Antennas 

Two  10-element  Yagi  antennas  are  used  (Taco,  Y-10JM)  which 
have  a  gain  of  10-12  db  above  that  of  an  isotropic  radiator. 
Although  the  nominal  impedance  of  these  antennas  is  50  ohms, 
it  was  found  that  at  221  me  the  terminal  impedance  was  ohms, 
necessitating  pi-matching  transformers  to  match  the  antennas  to 
the  30  ohm  transmission  line.  Each  antenna  is  mounted  to  a  polar 
axis  and  driven  by  a  sidereal  clock  drive.  RG8/U  coaxial  cables 
are  used  at  present  above  ground  between  the  antennas  and  the 
preamplifiers.  The  initial  antenna  separation  is  approximately 
109  feet. 

Pre- amplifiers 

The  pre-amplifiers  together  with  the  remainder  of  the  systems 
are  housed  in  a  24-foot  trailer  central  to  the  antenna  baseline. 
The  pre-amplifiers  are  of  conventional  design  consisting  of  two 
stages  of  grounded  grid  7077  planar  triodes,  followed  by  a  7077 
mixer  and  a  single  6AK5  IF  amplifier.  Using  an  "L"  type  (14) 
input  transformer  to  the  pre-amplifiers,  the  transfer  ratio  was 
experimentally  adjusted  to  obtain  an  optimum  noise  figure.  A 
noise  figure  in  the  range  of  4.5~5«0  db  can  be  maintained  in 
routine  operation  of  the  equipment. 

Allowing  for  cable  losses  and  an  antenna  gain  of  10  db,  the 
minimum  detectable  level  (using  the  criteria  given  by  Ityle  and 
Vonberg,  (15))  is  approximately  10  p  wm  ^cps"  for  a  receiver 
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bandwidth  of  600  me,  a  recording  time  constant  of  0.5  second, 
and  a  4.8  db  noise  figure.  Using  the  same  antennas  but  placing 
the  pre-amplifiers  at  the  antennas  would  result  in  an  improvement 
in  the  sensitivity  by  a  factor  of  2.5. 

I.  F,  amplifier 

A  standard  commercial  30  me  amplifier  manufactured  by  Lei  is 
used  in  the  system.  Over-all  IF  bandwidths  of  125  and  640  kc  are 
available  from  a  crystal  and  a  conventional  L/C  filter  respectively. 
Local  oscillator  system 

The  phase  lock  system  is  a  closed-loop  feedback  system  which 
is  used  to  produce  two  local  oscillator  signals  at  approximately 
251  me  which  are  separated  by  a  fixed  1  kc  difference  frequency. 

A  block  diagram  of  the  phase-lock  system  is  given  in  Figure  2.2. 

As  shown  in  the  figure,  the  two  fundamental  frequencies 
are  generated  in  fifth  overtone  crystal  controlled  Butler 
oscillator  circuits,  multiplied  in  doubler  circuits,  and  mixed 
to  obtain  the  difference  frequency.  The  1  kc  difference  signal 
is  amplified  and  fed  to  a  synchronous  detector  which  is  referenced 
to  a  1  kc  signal  obtained  from  a  reed  controlled  oscillator. 

The  oscillation  frequency  in  one  of  the  Butler  oscillators  is 
controlled  over  a  limited  range  by  the  dc  control  signal  formed 
at  the  output  of  the  synchronous  detector  as  shown  in  the  block 
diagram. 

A  circuit  schematic  of  the  voltage  controlled  crystal  oscil¬ 
lator  circuit  is  given  in  Figure  2.5.  As  shov.n,  non-linear 
capacitors  are  used  in  the  plate  circuit  of  V^.  A  change  in  the 
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BOEING  PHASE  LOCK  CONTROL  SYSTEM 


VOLTAGE  PULLED  CRYSTAL  CONTROLLED  BUTLER  OSCILLATOR  CIRCUIT 


dc  potential  applied  to  non-linear  capacitors  in  the  plate  circuit 
of  will  produce  a  change  in  capacity  and  result  in  a  phase 
shift  in  the  oscillator  feedback  loop.  The  phase  shift  results 
in  a  small  change  in  the  oscillator  frequency  to  achieve  phase- 
lock  equilibrium. 

The  polarity  of  the  control  signal  applied  to  the  capacitors 
is  such  that  nearly  a  phase  quadrature  condition  will  be  main¬ 
tained  (under  closed-loop  conditions)  between  the  two  input  signals 
at  the  synchronous  detector.  Since  the  1  kc  reference  signal  used 
in  the  lobe  sweep  system  is  obtained  from  mixing  the  two  local 
oscillator  signals,  it  should  be  clear  that  a  fixed  phase  relation¬ 
ship  is  maintained  between  this  reference  signal  and  that  resulting 
from  a  source  fixed  relative  to  the  antenna  pattern. 

Since  it  is  possible  to  achieve  phase  lock  at  odd  multiples 
of  the  1  kc  reed  oscillator  frequency,  some  initial  adjustment 
is  required  at  the  time  the  system  is  placed  in  operation.  It 
was  found  experimentally  that  the  control  range  was  sufficient 
to  maintain  phase  lock  without  the  use  of  temperature  stabilised 
ovens  for  the  crystals. 

The  frequency  of  the  local  oscillators  is  a  function  of  the 
inherent  stability  of  the  Butler  oscillators  while  the  difference 
frequency  stability  is  a  function  of  the  stability  of  the  reed 
controlled  oscillator  (nominally  ±1.5  cps).  The  difference 
frequency  will  remain  fixed  independent  of  frequency  drifts  of 
the  crystal  oscillator  frequencies  as  long  as  the  drift  is  within 
the  pull-in  range  of  the  voltage  controlled  oscillator. 
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The  use  of  a  crystal  controlled  oscillator  as  the  basic  frequency 
source  takes  advantage  of  the  inherently  high  "Q"  of  the  crystal, 
thereby  minimizing  problems  associated  with  spurious  oscillations, 
multiple  oscillation  modes  and  residual  frequency  modulation 
components.  A  high  degree  of  shielding  was  required  to  provide 
adequate  isolation  between  the  two  local  oscillator  channels. 

Cross  coupling  between  the  two  channels  will  result  in  a  residual 
1  kc  component  in  the  system  and  can  also  result  in  spurious  fre¬ 
quency  modulation  of  the  1  kc  difference  frequency. 

AGC  system  and  intensity  recorders 

A  block  diagram  outlining  the  AGC  system  is  shown  in  Figure  2.4A. 
Referring  to  the  figure,  it  is  seen  that  the  1  kc  modulated  noise 
signal  is  amplified  in  the  main  IF  amplifier,  square  law  detected 
and  passed  through  two  filters  30  and  3  cps  wide.  The  first  filter 
is  a  conventional  passive  L/C  filter.  The  second  filter  is 
approximately  3  cps  wide  and  employs  a  bridged-"T"  network  in  a 
feedback  amplifier  (l6). 

As  shown  in  the  figure,  the  1  kc  component  of  the  signal  after 
filtering  is  threshold  detected,  chopped  at  a  417  cps  rate,  amplified, 
and  synchronously  detected  and  filtered  to  obtain  the  dc  gain  control 
signal  which  is  fed  back  to  the  main  IF  amplifier.  The  use  of  the 
chopper  and  synchronous  detector  allows  for  versatility  and  ensures 
high  dc  stability  of  the  AGC  delay  level.  The  gains  of  the  various 
amplifiers  are  set  so  that  the  low  range  recorder  indicates  approxi¬ 
mately  20  percent  full  scale  for  quiet  sun  levels  and  has  a  uniform 
logarithmic  gradient  from  this  point  to  80  percent  full  scale,  a  level 
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FLIP-HOP 


10  db  above  the  quiet  sun.  The  deflection  of  the  extended  range 
recorder  shown  in  the  block  diagram  is  a  logarithmic  function  of 
input  coherent  noise  level  above  the  AGC  delay.  Full  scale 
deflection  of  this  recorder  corresponds  to  increases  of  approxi¬ 
mately  35  db  above  that  of  the  quiet  sun. 

1  kc  phase  meter 

The  1  kc  phase  meter  converts  the  phase  difference  between 
the  1  kc  signal  and  reference  to  an  analog  voltage  output.  The 
scheme  used  to  obtain  an  analog  measurement  of  phase  difference 
has  been  described  before  (2S)  and  for  the  present  circuit 
consists  of  triggering  a  flip-flop  to  turn  on  and  off  at  the 
time  of  the  negatively-going  zero  cross-over  for  the  1  kc  reference 
and  signal  respectively.  The  integrated  output  of  the  flip-flop 
using  this  scheme  varies  linearly  with  phase  difference,  repeating 
itself  every  2*  radians.  An  input  phase  variation  at  a  constant 
rate  will  therefore  produce  a  sawtooth  waveform  at  the  output  of 
the  phase  meter. 

The  flip-flop  is  shown  in  the  block  diagram  of  the  phase  meter 
in  Figure  2.4B.  The  two  trigger  pulses  are  derived  from  the  1  kc 
reference  signal  and  the  1  kc  coherent  noise  signal  which  are 
amplified  and  limited,  using  cathode-coupled  double-ended  clipper 
circuits  (17).  The  1  kc  clipped  signal  output  from  each  channel 
is  fed  to  bi-static  transistor  Schmitt  triggers  which  fire  at 
the  zero  crossing  of  the  input  signals.  The  square  wave  outputs 
from  the  Schmitt  circuits  are  differentiated  to  form  the  triggers 
for  the  flip-flop. 
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Using  the  above  scheme  it  is  found  that  the  output  of  the 
analog  phase  recording  is  comparatively  insensitive  to  changes 
in  input  power  level.  The  output  deflection  changes  a  maximum 
of  1  percent  over  an  input  power  range  of  0  -  35  db  increase 
above  that  of  the  quiet  sun. 

The  signal  input  to  the  phase  meter,  as  shown  in  Figure  2.1, 
comes  from  a  1  kc  filter,  which  has  a  bandwidth  of  30  cps.  The 
30  cps  bandwidth  is  sufficiently  narrow  to  prevent  multiple 
zero  crossovers  in  the  signal.  The  presence  of  internal  noise 
will  result  in  phase  modulation  only  of  the  signal  at  the  output 
of  the  limiter  shown  in  the  phase  meter  block  diagram.  The 
phase  modulation  will  result  in  statistical  fluctuations  of  the 
analog  phase  output,  with  the  mean  output  representing  the  true 
phase  difference  between  the  reference  and  1  kc  signal  in  the 
absence  of  noise. 

Calibration  system 

The  system  is  calibrated  to  measure  power  levels  relative  to 
an  arbitrary  '  evel  output  from  a  5722  saturated  noise  diode.  This 
calibration  le  '1  is  chosen  to  represent  a  level  of  power  output 
close  to  that  expected  for  the  quiet  sun.  3y  comparing  this  power 
level  with  the  power  received  from  Cassiopeia,  the  system  is 
calibrated  to  measure  flux  density  in  units  of  watts  per  square 
meter  per  cycle  per  second.  From  (18),  the  level  received  from 
Cassiopeia  was  taken  to  be  1.0  x  10  wm  ^(cps)  at  a  frequency 
of  221  me.  The  zero  db  or  "quiet  sun"  level  is  ten  times  the 
level  from  Cassiopeia. 
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Two  levels  of  power  are  inserted  at  the  input  to  the  preamplifiers 
each  hour.  These  levels  correspond  to  that  of  the  quiet  sun  and  the 
quiet  sun  plus  10  db.  The  output  of  a  wide-band  220  me  amplifier 
serves  to  provide  the  +10  db  calibration  level.  A  daily  check  is 
made  to  determine  the  variations  of  the  noise  output  of  the  ampli¬ 
fier  relative  to  the  noise  calibration  diode.  Stability 
tests  and  spot  daily  checks  have  established  that  both  the  long 
and  short  term  output  noise  stability  from  the  amplifier  are 
within  1  db.  The  amplifier  has  approximately  40  db  gain  and 
consists  of  4  stages  of  grounded  grid  7077' s.  This  comparatively 
large  available  output  noise  power  is  used  in  checking  the 
response  and  calibration  of  the  system  over  a  35  db  range. 

Calibration  levels  of  20,  30,  and  35  db  are  inserted  into  the 
system  once  per  day.  Stability  tests  are  performed  on  a  daily 
basis. 

All  the  rf  switching  is  accomplished  through  the  use  of 
motor  driven  coaxial  switches.  A  variable  line  stretcher  is 
used  to  cancel  out  the  effect  of  the  cross  coupled  components, 
via  the  common  noise  calibration  source,  between  the  two  pre¬ 
amplifiers.  In  practice,  the  variable  line  stretcher  is  ex¬ 
perimentally  adjusted  for  a  null  output  on  the  intensity  recording. 

The  null  is  sufficiently  below  the  quiet  sun  level  in  the  system 
that  no  appreciable  error  is  introduced  in  intensity  measurements. 
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PART  III 


CONSIDERATIONS  RELATING  TO  POSITION  MEASUREMENTS 

Introduction 

This  section  of  the  report  deals  with  the  geometrical 
considerations  relevant  to  the  determination  of  burst  position. 

The  system  of  coordinates  is  so  chosen  that  the  results  can  be 
applied  to  an  interferometer  of  arbitrary  orientation.  Hie  derived 
results  are  also  used  in  the  analysis  of  the  phase  compensator  in 
the  next  part  of  this  report.  Hie  received  data  are  discussed 
and  it  is  shown  how  they  are  interpreted  to  provide  burst  position 
measurements.  The  error  in  position  measurements  resulting  from 
the  fact  that  the  radio  and  optical  center  of  the  sun  may  not 
coincide  is  investigated. 

Solar  coordinate  system 

With  respect  to  a  fixed  observer  a  point  on  the  celestial 
sphere  can  be  specified  by  its  declination  (6)  and  its  hour 
angle  (HA)  (Figure  3.1).  By  virtue  of  the  rotation  of  the  earth 
the  hour  angle  increases  continuously.  The  declination  of 
celestial  objects  is  unaffected  by  the  earth's  rotation;  however, 
the  declination  of  the  sun  changes  slowly  due  to  the  revolution 
of  the  earth  around  that  body. 

It  will  be  convenient  to  specify  a  position  on  the  sun  with 
respect  to  some  reference  point  called  the  "center."  Hie  axes 
for  the  reference  system  will  be  the  hour  circle  and  the  parallel 
of  declination  intersecting  perpendicularly  at  this  center. 
Distances  along  the  hour  circle  will  be  indicated  by  A 6.  However, 
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POUR  AXIS 


SPECIFICATION  OF  A  POINT  ON  THE  CELESTIAL 
SPHERE  RY  ITS  DECLINATION  ANGLE  AND  HOUR  ANGLE  . 

FIG.  3.1 
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since  the  increment  of  hour  angle  subtended  by  the  solar  diameter 

increases  with  declination,  distances  along  the  parallel  of 

declination  will  be  indicated  by  AHA  cos  6q  so  that  angular 

distances  will  be  the  same  in  both  directions.  The  declination 

6  is  that  of  the  sun's  center, 
o 

In  general,  the  locus  of  possible  source  positions  from 
interferometer  data  is  a  straight  line,  over  angular  distances 
of  the  order  of  the  sun's  diameter.  This  line  can  conveniently 
be  specified  by  the  length  r  and  polar  angle  t  of  its  normal 
vector  from  the  center  (Figure  3.2).  tis  measured  counterclockwise 
from  the  west. 


SOUR  COORDINATES 


FIG.  3.2 
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Baseline  coordinates 


It  is  desired  to  calculate  the  position  in  solar  coordinates 
of  activity  centers  on  the  sun  from  the  radio  data.  This  cal¬ 
culation  is  dependent  on  the  length  and  orientation  of  the 
interferometer  baseline,  which  must  be  specified. 

Given  two  antennas  A  and  B,  with  inputs  X  cos  (tot  +  0^)  and 
Y  cos  (cot  +  0  )  respectively,  it  can  be  shown  that  the  received 

y 

signal  input  to  the  phase  meter  is  proportional  to 

XY  cos  (to  t  —  0  +  0  +  0  ) 

a  x  y  a 

where  0a  is  a  constant  relative  phase  difference  between  channels, 
and  '•)  is  the  lobe  sweep  frequency.  The  antenna  A  is  defined  as 
the  antenna  whose  phase  subtracts,  while  3  is  the  antenna  whose 
phase  adds, so  that  the  phase  meter  reads  directly  the  quantity 
(-0X  +  0y  +  0a). 

Define  the  baseline  as  a  vector  extending  from  A  to  B.  Then 
this  vector  is  the  parameter  necessary  for  conversion  from  radio 
data  to  solar  coordinates. 

The  baseline  vector  will  be  specified  by  its  length,  its 
inclination  y  to  the  polar  axis,  and  its  position  angle  {  measured 
from  the  north  point  in  the  equatorial  plane  westward  to  the 
projection  of  the  vector  in  the  equatorial  plane  (Figure  3.5). 
Derivation  of  the  incident  angle 

The  phase  difference  between  the  two  interferometer  antennas 
depends  on  the  incident  angle  9  between  the  plane  normal  to  the 
baseline  and  the  incident  ray.  This  angle  will  be  considered 
positive  if  the  ray  path  to  the  B  antenna  is  shorter  than  the 
path  to  the  A  antenna,  and  negative  otherwise. 
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POLAR  AXIS 


SPECIFICATION  OF  IASELINE  ORIENTATION 
FIG.  3.3 
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Figure  3.4  is  a  geometrical  diagram  of  the  relationships 
involved  between  sun  position  and  incident  angle*  This  is  a 
diagram  formed  by  looking  approximately  normally  at  both  the 
polar  axis  and  the  baseline.  From  this  diagram  it  can  be  seen 
that 


3.1  sin  6  =  sin  6  cos  y  -  cos  6  cos  (HA  +  §)  sin  y 

From  the  angle  0,  the  quantity  0  -  0  can  be  found  (Appendix 

y  ^ 

A).  The  ray  path  difference  between  the  antennas  is  d  sin  0. 

2  ^  j 

This  represents  a  phase  difference  of  sin  9. 

Therefore,  0  -  0  =  sin  9 

y  x  * 

Substituting  (3.1)  gives 

3.2  0  -  0  *  [sin  6  cos  y  -  cos  6  cos  (HA  +5)  sin  y3 

y  ^ 

For  a  horixontal  east-west  baseline  r  *  5  *  Thus  from 
(3.2)  we  get 

3.3  0  -  0  c  ain  HA  cos  6 

y  x 

for  an  east-west  baseline. 

Conversion  to  solar  coordinates 

It  is  desired  to  convert  a  change  in  0  -  0  to  the  solar 

jr  * 

coordinates  As  and  AHA  cos  6  •  Partially  differentiating  (3*2) 

o 

gives 

A(0y~0x)»^jffi  (cos6ocosY)A6+(sin6ocos(HA+5)sinY)A6+(cos6osin(HA+5)sinY)AHA] 

Let  A  HA  cos  6  =  AH.  Then 

o 
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A(0  -  0  )*  A0  s  [(cos  6  cos  y  +  ain  6  coa  (HA  +  ?)  ain  y)  A6  + 

y  x  a  o 

ain  (HA  +  ?)  ain  y  Ah] 

Thus,  for  a  given  phase  meter  change  there  ia  obtained  the  equation 
of  a  atraight  line  specifying  the  location  of  the  source.  The 
equation  ia  of  the  form 


A0  =  AA6  +  BAH 


where 


3,4  A  s  ^i^Ccos  6  coa  y  +  ain  6  coa  (HA  +  ?)  ain  y] 

AO  O 


3*5  B  e  — sin  (HA  +5)  sin  y 

The  intercepts  on  the  AH  and  A6  axes  are  ^  and  ^respectively 
(see  Figure  3.5). 

From  Figure  3.5  can  be  found  the  parameters  r  and  *  in  terms  of 
A,  B,  end  A0: 


3.6  ain  t  «  tS 


w* 


3.7  coa  t  = 


Therefore 


3.8  t  -  tan"1  | 

The  angle  is  still  ambiguous  by  l80°  aa  so  defined.  To  resolve 
this  ambiguity,  the  sign  of  the  y  intercept,  can  be  employed. 
The  following  rule  will  be  used: 

Rule:  Assume  $  in  the  first  or  second  quadrant.  Compute 
If  this  ia  negative,  add  ft  to  the  angle. 
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AS 


LOCUS  IN  SOLAR  COORDINATES 
FIG.  3.5 
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Squaring  (3.6)  and  (3.7)  and  adding  gives: 


1 

(A0)2 


[f|j  +  G$] 

r2  (A2  +  B2),  or, 


3.9  r  .  ■ 

7  A2  +  B2 

Datarnination  of  burat  poaitione  from  the  racordad  measurements 

As  was  discussed  in  Section  II  tha  received  data  consist  of 
two  principal  recordings  which  are: 

1 —  An  analog  presentation  of  the  phase  difference  of  the 
received  signal  between  the  two  receiving  antennas. 

2—  An  output  which  is  a  measure  of  the  intensity  of  the 
total  available  coherent  noise  power. 

With  square-law  detection  the  recorded  phase  and  amplitude 
are  the  resultant  of  the  superposition  of  all  coherent  noise 
sources  within  the  primary  beam  of  the  receiving  antennas.  Thus 
the  quiet  sun  and  an  enhanced  component  of  emission  will  both 
produce  a  component  of  phase  and  amplitude  on  the  two  respective 
recordings.  The  phase  recording  will  be  displaced  from  the  level 
due  to  the  quiet  sun  if  the  enhanced  component  is  not  located  at 
the  position  of  the  radio  center  of  gravity. 

A  typical  recording  is  shown  in  Figure  3.6.  The  top  trace 
is  the  intensity  recording.  The  bottom  record  is  an  analog 
recording  of  the  phase  measurement.  The  saw-tooth  appearance 
of  the  phase  recording  results  from  the  sun  moving  across  the 
interferometer  lobe  pattern.  Each  cycle  corresponds  to  a  change 
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in  path  difference  of  the  received  signal  between  the  antennas  of 
one  wave  length.  All  of  the  measurements  to  date  have  been 
recorded  in  the  manner  shown  in  Figure  3*6.  The  phase  compensator, 
which  will  be  incorporated  into  the  system  shortly,  will  have  the 
function  of  inserting  into  the  system  a  variable  phase  nearly 
equal  to  that  predicted  for  the  optical  sun.  This  will  ideally 
convert  the  saw-tooth  phase  recording  into  a  straight  line. 

At  present  the  measurements  of  burst  positions  of  enhanced 
emissions  are  made  relative  to  the  radio  center  of  gravity, 
which,  in  general,  does  not  coincide  with  the  optical  center 
of  the  sun. 

The  resultant  amplitude  and  phase  of  the  signal  into  the 
phase  meter  during  the  burst  will  be  the  superposition  of  two 
independent  sine  waves  resulting  from  that  of  the  pre-burst 
radio  sun  and  that  from  the  burst.  To  determine  the  position 
of  the  burst  relative  to  the  radio  center  of  gravity,  it  is 
required  to  determine  the  electrical  phase  angle  between  the 
signal  arising  from  the  burst  and  that  of  the  signal  resulting 
from  the  sun  in  the  absence  of  the  burst.  The  relationship 
between  the  measured  phase  shift,  due  to  the  burst  plus  the 
pre-burst  sun,  and  the  phase  shift  of  the  burst  component  is 
derived  in  Appendix  A  and  is 

P2  ,  , 

jr  8in  Ax 

3.10  L0  »  tan  -1  - - 

*2 

cos  AX  -1 
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where  Pg/P^  *  measured  ratio  of  intensities  before  and  after 
burst 

Ax  =  measured  phase  deviation  resulting  from  burst 
plus  pre-burst  radio  sun 
A0  =  relative  phase  of  the  burst 
Referring  again  to  Figure  3.6,  showing  the  analog  phase  output, 
it  is  seen  that  Ax  is  obtained  simply  by  noting  the  change  in  phase 
deflection  at  the  time  of  the  burst.  The  phase  shift  Ax  is  the 
phase  shift  to  the  new  radio  center  of  gravity,  somewhere  between 
the  pre-burst  center  of  gravity  and  the  burst.  The  phase  shift 
A0,  which  cannot  be  measured,  is  the  phase  shift  due  to  the 
burst  alone. 

Combining  (3.9)  and  (3*10),  the  resultant  burst  position  for 
an  arbitrary  baseline  orientation  is 


3.11 


r 


(Pg/^)  8in  Ax 

00.  ix  -  1 

UTW71 


For  a  horizontal  east-west  baseline  r  ■  5  ■  */2  (see  Figure  3.3) 
and  from  (3.4)  and  (3.3)  we  have 


A  a  sin  sin  HA 

B  a  — cos  HA 


Therefore  from  (3.9)  and  (3.10)  for  a  horizontal  east-west 
baseline,  the  burst  position  is  given  by 


3.12 


r 


\ 

2nd 


1  (P2/P1)  8±n  AX 

t,n~  tyiyr  ~  -  i 

(1  -  .in2  HA  CO.2  6  )1/Z 

o 
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Similarly,  the  orientation  angle  for  a  horizontal  east-west 


baseline  from  (3*8)  is 

3.13  ♦  ■  tan"^  (-sin  6  tan  HA) 

o 


To  illustrate  the  method  used  in  the  determination  of  burst 
positions,  consider  the  7.4  db  isolated  burst  shown  in  the  received 
data  given  in  Figure  3.8.  The  relevant  data  relating  to  the  burst 
are  as  follows: 

declination  of  the  sun  ■  6  ■  21°  18 ' 

hour  angle  of  the  sun  at  the  time  of  the  burst  «  8.25° 

p2/pi  «  7.4  db  «  5.5 

measured  phase  shift  «  4x  *  29° 

operating  frequency  «  221.34  me 

antenna  separation  (horizontal  east-west)  »  33.2  meters  ■ 
109.0  feet 

From  equation  (3,10)  the  electrical  phase  shift  of  the  burst  is 
calculated  to  be  34.93°.  Using  equation  (3.12)  tire  calculated 
burst  position  is  13.6' •  From  equation  (3.13)  the  orientation 
angle  is  computed  to  be  -2°  34'.  The  burst  is  thus  located  to  the 
west  of  the  radio  center  of  gravity. 

Error  in  the  determination  of  burst  position  resulting  from  an 
error  in  the  assumed  position  of  the  radio  center  of  the  sun 
From  (3.10)  and  (3.12),  it  is  seen  that  the  burst  position 
given  for  a  horizontal  east-west  baseline  is 


3.14 


r 


\  _  AflT 

^  (1  -  sin2  HA  cos2  6  )X/^ 

o 
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and,  from  (3*1)  this  becomes 


3-15  r "  m 

The  incident  angle  0  ia  the  angle  of  the  pre-burst  radio  center 
of  gravity  of  the  sun  which  is  assumed  to  ooineide  with  the  solar 
optical  center*  An  error  in  the  incident  angle  of  the  radio  center 
of  gravity  will  result  in  an  error  in  the  computed  burst  position. 
From  (3*13),  the  fractional  error  is 

to  eo‘  *.  . 

r  "  cos  (0  +  A©) 

0 

From  a  Taylor  expansion  of  the  denominator  for  AO  small  we 

have 

5-16 

The  error  turns  out  to  be  negligible  except  for  large  incident 
angles.  Given  an  incident  angle  of  83  degrees  and  an  error  in  the 
incident  angle  of  13' ,  the  resulting  error  is  3.2  percent.  Obser¬ 
vations  reported  by  Boischot  (10)  at  169  me  indicate  that  the 
maximum  deviation  of  the  center  of  gravity  from  the  optical 
center  will  not  exceed  1  solar  radius. 
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PART  IV 


THE  PHASE  COMPENSATOR 


Introduction 

Since  the  phase  compensator  represents  a  major  effort  in  the 
solar  research  program  and  is  an  innovation  in  interferometer 
instrumentation,  a  section  of  this  report  is  devoted  to  the  dis¬ 
cussion  of  this  new  device.  The  compensator  is  currently  under¬ 
going  tests  and  has  not  as  yet  been  incorporated  into  the  system 
for  routine  solar  observations. 

ji s  was  mentioned  previously,  the  phase  compensator  will 
improve  the  performance  of  the  system  in  two  respects:  it  will 
remove  the  phase  discontinuity  from  the  present  sawtooth  waveform 
output  of  the  phasemeter,  and  it  will  enable  the  phasemeter  to  indicate 
the  variation  in  position  of  the  radio  center  of  gravity  of  the  sun. 

The  phase  compensator  will  accomplish  these  by  generating  and 
inserting  into  the  system  a  time-varying  phase  which  is  nearly 
equal  to  that  predicted  for  the  optical  center  of  the  sun  as  it 
moves  relative  to  the  antenna  baseline.  This  generated  time- 
varying  phase  will  be  subtracted  from  the  time-varying  phase  of 
the  received  signal  so  that  ideally  the  present  sawtooth  waveform 
of  the  phasemeter  will  be  converted  to  a  straight  line. 

This  section  of  the  report  discusses  the  principle  of  operations 
of  the  phase  compensator  and  the  function  of  its  various  elements. 

The  last  part  of  this  section  considers  the  error  which  results 
from  the  fact  that  the  compensator  does  not  correct  for  changes 
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in  the  declination  of  the  sun,  or  in  the  equation  of  time.  The 
error  is  analyzed  in  the  special  case  of  an  east-west  interferometer 
baseline.  The  calibration  procedure  for  the  compensator  is  also 
discussed. 

Preliminary  considerations 

It  will  be  assumed  for  simplicity  that  the  optical  and  radio 
centers  of  the  sun  coincide.  Equation  (3.2)  can  then  be  written 

4.1  0  -  0  b  £i£[sin  6  cos  y  -  cos  6  cos  (H.-*  +  §)  sin  y] 

y  x  a  o  o  o 

where  6  and  HA  are  the  declination  and  the  hour  angle  of  the 
o  o 

optical  center  of  the  sun,  respectively.  The  phase  meter  reads 
this  value  directly. 

Under  normal  operation  of  the  phasemeter,  the  largest  single 
contribution  to  the  recorder  output  is  due  to  the  second  term  of 
(4.1).  This  completely  masks  any  long  term  variations  of  the 
sun's  radio  center,  and  in  addition  makes  more  difficult  the 
measurement  of  short  term  variations.  It  is  desirable,  therefore, 
to  subtract  this  term  out.  However,  this  must  be  done  in  such  a 
manner  that  the  process  of  subtraction  does  not  in  itself  introduce 
appreciable  error.  This  rules  out  the  use  of  an  analog  device, 
which  would  permit  the  accumulation  of  error  with  time. 

If  the  quantity  L-  -j--  cos  6q  cos  (HAq  +5)  sin  y]  is  subtracted 
from  the  phase  meter  reading,  then  the  phase  meter  will  read 

4.2  (0  -  0  )'  ■  sin  6  cos  y 

y  x  a  o 

which  is  almost  a  straight  line  as  a  function  of  time.  Any  observed 
deviations  from  a  straight  line  would  be  due  to  changes  in  the 
position  of  the  source  alone. 
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The  quantity  [-  cos  6  sin  y  cos  (HA  +  ?)]  is  a  quasi- 

n  O  O 

sine  wave,  with  amplitude  [-  cos  &q  sin  yl  *  period  of  one  day, 
and  phase  5.  The  term  "quasi"  is  used  because  cos  &q  varies 
slightly  over  the  course  of  a  day  and  HA  does  not  vary  linearly 
with  time  because  of  the  equation  of  time  perturbation. 

Summary  of  the  operating  principles  of  the  phase  compensator 

The  phase  compensator  inserts  into  the  reference  signal  input 
to  the  phase  meter  a  phase  shift  which  varies  as  a  sine  wave. 

This  is  accomplished  by  storing  on  paper  tape  the  differences 
of  sin  HA,  for  15"  increments  of  arc,  over  the  range  of  H.t  from 
0°  to  90°.  Each  second  of  time,  a  paper  tape  reader  reads  out 
this  incremental  difference.  The  difference  is  converted  to  a 
pulse  train  in  which  the  number  of  pulses  is  proportional  to  the 
difference  of  sin  HA  read  out.  The  resulting  pulse  train  is 
fed  to  a  ratio  counter  which  takes  out  a  preset  fraction  of 
the  pulses  sent  to  it.  This  is  the  method  for  setting  the 
amplitude  of  the  sine  wave.  The  output  pulses  from  the  fractional 
divider  are  used  to  drive  a  stepping  motor  which  in  turn  drives 
a  resolver,  where  the  subtraction  mentioned  earlier  takes  place. 

The  tape  reader  is  synchronized  to  mean  solar  time  with  the 
origin  of  the  sine  wave  set  so  that  the  time  of  maximum  rate  of 
change  corresponds  to  meridian  passage  of  the  sun  (for  an  east- 
west  baseline).  The  operation  is  continuous  with  the  tape 
reversing  direction  each  90  degrees  of  hour  angle. 

For  a  radio  star  source  and  a  time  standard  synchronized  to 
sidereal  time  instead  of  mean  solar  time,  no  day-to-day  correction 
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would  be  required,  as  the  inserted  phase  variation  would  exactly 
duplicate  that  caused  by  the  rotation  of  the  earth*  Since  the 
equation  of  time  and  declination  of  the  sun  change  over  a  24-hour 
period  and  the  amplitude  of  the  sine  phase  curve  inserted  by  the 
compensator  is  constant  over  any  one  24-hour  period,  it  is  clear 
that  an  error  will  develop. 

The  phase  compensator  equipment 

The  phase  compensator  consists  of  nine  parts:  (l)  a  paper 
tape  input,  on  which  is  punched  the  sine  wave  data,  (2)  a  timing 
system,  which  controls  the  pulse  flow  through  the  system,  (3)  a 
data  counter,  which  converts  the  tape  data  to  a  pulse  sequence, 

(4)  a  ratio  counter,  which  provides  a  fine  control  of  the  amplitude 
of  the  sine  wave  by  selecting  a  pre-determined  fraction  of  the 
pulses  in  the  sequence  for  removal,  (5)  an  inhibit  circuit,  which 
inhibits  the  pre-determined  fraction  of  input  pulses  on  command 
from  the  ratio  counter,  (6)  a  divider  chain,  which  provides 
coarse  control  of  the  sine-wave  amplitude,  (7)  a  stepping  motor, 
which  converts  the  pulse  sequence  to  rotary  motion,  (8)  a  gear 
step-down  box,  and  (9)  a  resolver,  for  converting  the  rotary  motion 
to  a  phase  subtraction.  These  elements  will  be  discussed  in  turn. 
Figure  4.1  is  a  block  diagram  of  the  system. 

Tape  input 

On  the  tape  are  punched  values  of  the  fifth  decimal  place 
differences  between  values  of  a  sine  function  for  consecutive  15" 
intervals,  that  is,  the  incremental  difference  punched  at  a  given 
position  is  the  difference  between  the  value  of  the  sine  at  the 
given  position  and  the  sum  of  the  differences  already  punched. 
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HOCK  DIAGRAM  OF  PHASE  COMPENSATOR 


The  fifth  place  differences  change  at  most  by  eight  digits  for 
an  increment  of  15".  Thus,  sin  0°  ■  0,  sin  15"  ■  0.000074.  The 
tape  is  advanced  one  position  each  second  of  mean  solar  time,  by 
a  pulse  from  the  timing  system.  At  each  position  there  is  read 
out  from  the  tape  the  sine  difference  value  punched  there,  a 
number  from  zero  to  eight.  This  number  is  then  sent  to  the  data 
counter. 

Timing  system 

The  timing  system  consists  of  a  60  cycle  per  second  time 
standard  and  associated  pulse-forming  circuitry.  The  timing 
system  provides  a  one-second  synchronizing  pulse  which  reads 
out  and  advances  the  tape  once  a  second.  This  is  accomplished 
by  a  six  stage  binary  counter  which,  once  a  cycle,  by  means  of 
switching  circuitry,  has  a  count  of  4  added  to  it.  Thus  the 
complete  cycle  is  60  counts  instead  of  64. 

The  timing  system  also  provides  three  pulse  sequences  at 
the  50  cps  repetition  rate,  synchronized  with  the  one-second 
pulse,  but  with  three  different  delays,  achieved  with  HC  delay 
networks  and  clamping  diodes.  These  control  the  switching  in 
various  parts  of  the  compensator. 

7 

The  time  standard  is  accurate  to  5  parts  in  10  per  month. 

It  is  unaffected  by  line  frequency  variations. 

Data  counter 

The  data  counter  is  the  heart  of  the  sine  generation  process. 
The  one-second  synchronizing  pulse  passes  through  the  holes  in 
the  paper  tape  and  reads  a  number  into  the  data  counter.  The 
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gating  is  such  that  the  number  read  into  the  counter  is  -7  times 
the  number  in  the  tape.  Thus,  a  2  becomes  a  -14.  As  soon  as  the 
data  counter  switching  circuitry  senses  that  its  contents  are  not 
zero,  the  data  counter  begins  to  count  at  a  60  cps  rate,  meanwhile 
providing  an  output  pulse  for  every  input  pulse  from  the  timing 
system.  When  the  contents  of  the  data  counter  are  again  zero, 
i. e. ,  when  the  number  of  input  pulses  from  the  timing  system 
equals  the  number  received  from  the  tape,  it  stops  and  awaits  the 
next  input  from  the  tape. 

Thus,  the  data  counter  generates  a  train  of  pulses  at  a  60  cps 
rate  equal  in  number  to  7  times  the  number  stored  on  the  tape.  At 
maximum  rate  of  change  of  the  sine  wave,  an  S  will  appear  on  the 
tape,  and  during  a  60-cycle  second,  the  data  counter  will  provide 
56  pulses  on  the  first  56  cycles. 

Ratio  counter 

The  ratio  counter  is  a  14  stage  binary  counter  which  selects 
a  pre-determined  fraction  of  pulses  from  the  data  counter  for 
removal.  The  first  stage  turns  "on"  for  every  other  input  pulse, 
the  second  stage  turns  on  for  every  fourth  input  pulse,  and  so 
forth.  Thus,  the  first  stage  can  be  gated  to  the  inhibit  circuit 
to  remove  one-half  the  input  pulses,  the  second  stage  to  remove 
one-fourth  the  input  pulses,  and  so  forth.  This  gating  is 
accomplished  with  the  use  of  fourteen  switches  which  are  set 
at  calibration  time.  Any  fraction  F  of  the  form 


F 


s  0  or  1 
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can  be  built  up  by  setting  the  switches  according  to  the  o^'s. 

If  the  j  switch  is  in  the  off  position,  this  switch  will  provide 
for  the  rejnoval  of  of  the  input  pulses. 

Under  consideration  is  the  use  of  a  ratio  device  such  &s  this 
for  a  solar  to  sidereal  time  converter.  By  using  enough  stages, 
any  desired  accuracy  could  be  obtained,  such  as  one  second  per 
year. 

Inhibit  circuit 

The  inhibit  circuit  is  where  the  actual  pulse  removed  takes 

place.  '.Then  the  position  in  the  ratio  counter  goes  on,  and 
til 

the  j  switch  is  off,  a  pulse  from  the  ratio  counter  sets  the 
inhibit  trigger.  After  a  delay,  the  pulse  arrives  from  the  data 
counter  and  is  inhibited.  Then  after  another  delay,  a  pulse 
arrives  to  reset  the  inhibit  trigger,  which  now  awaits  the  next 
command  from  the  ratio  counter. 

Divider  chain 

This  is  simply  a  seven-stage  binary  counter  which  reduces  the 
pulse  rate  by  — where  n  can  be  selected  by  means  of  a  rotary 
switch,  from  zero  to  seven.  The  number  n  is  fixed  for  a  given 
baseline  such  a  value  that,  for  a  given  total  fraction  of 

pulses  removed,  a  minimum  fraction  is  removed  in  the  ratio  counter 
and  a  maximum  fraction  in  the  divider  chain. 

3tep  motor,  gear  reduction,  resolver 

These  are  self-explanatory.  The  stepping  motor  rotates  18° 
per  input  pulse.  This  is  stepped  down  by  36:1  to  0,5°  per  input 
pulse  by  the  gear  train,  which  in  turn  rotates  the  resolver. 
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Setting  the  compensator  for  a  particular  set  of  parameters 

This  is  best  illustrated  with  an  example.  The  amplitude 

of  the  reference  sine  wave  which  is  to  be  subtracted  from  the 

2  rtd 

phase  sine  wave  of  the  received  signal  is  cos  6  sin  y,  in 
units  of  radians.  Assume  the  following  parameters: 


2nd 


T 


189.46 


6  *  cos"1  (0.92)  «  23° 


Y  =  5  «  */2 


Therefore  the  amplitude  is 

cos  6  sin  Y  *  174.3  radians  =  10,000  degrees. 


This  means  that  during  the  course  of  six  hours  from  mean  noon 
during  which  the  hour  angle  has  increased  by  90°  and  the  phase 
sine  wave  has  gone  from  aero  to  full  amplitude,  the  motor  must 
receive  20,000  pulses,  since  it  advances  0.3°  per  pulse.  The 

data  counter,  on  the  other  hand,  during  this  same  period  generates 

5  2 

a  sequence  of  pulses  equal  to  7  x  10  .  Thus,  a  reduction  of  ^ 
is  called  for.  This  reduction  can  be  written 


* 


The  reduction  of  called  the  coarse adjustment,  is  set  by 
tapping  off  the  output  of  the  fifth  stage  of  the  divider  chain. 
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The  ratio  of  ||  ■  .914286  is  Bet  in  the  ratio  counter.  The 
ratio  counter  can  be  set  to  inhibit  any  number  of  input  pulses 
during  the  course  of  its  21Zf  «  16,  384  count  cycle.  The  exact 
number  we  would  like  to  pass  is 

||  x  16,384  -  14,979.7 

But,  since  we  can  pass  only  an  integral  number  of  pulses,  we  must 
be  content  with  14,980.  Thus,  the  true  ratio  will  be 

lH§5  “  • 914307 

The  error  is  ~  .00002 

To  set  the  ratio  counter  to  pass  the  14,980  pulses,  first  write 
the  number  less  one  in  binary  form  as  follows: 

14,979  »  11101010000011 

Ihen,  if  the  binary  digits  are  thought  of  as  stages  of  the  ratio 
counter,  all  stages  with  a  1  should  be  in  the  on  (not  inhibit) 
position.  The  maximum  error  is 

_ _ _ 

Number  of  counts  passed 

The  greater  the  number  of  counts,  the  smaller  the  percent  error. 

The  divider  chain  permits  the  passing  of  greater  than  ^  x  16,384 
pulses.  Thus,  the  greatest  error  is  of  the  orde’*  of  .00006. 

Phase  compensator  error 

The  phase  of  the  1  kc  signal  input  to  the  phase  meter  is  equal 
to  the  phase  difference  between  the  receiving  antennas. 
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From  (4*1)  the  phase  of  tho  1  ko  signal  for  an  east-west  baaalina 
(where  r  and  g  ■  */2)  is  given  by 

4.3  *slg  ■  TT  •*>  BAo  «• s. 

Taking  into  aeoonnt  a  ehangs  in  daelination  aqual  to  A6q  over 
a  24-hour  period,  and  with  tho  tins  t  given  in  hours,  (4*3)  bocoass 

4.4  ^sig  *  coa  (5i +  "2f"}  Bla  HA0 

whoro  6*  is  tho  doclination  of  tho  oontor  of  tho  sun  at  t  ■  o. 
o 

Prior  to  calibration,  tho  compensator  insorts  a  nearly  identical 
time-varying  phase  into  tho  roforoneo  signal  input  to  tho  phase 
motor.  Since  tho  amplitude  of  tho  phase  sine  wave  output  from 
tho  eospensator  can  bo  varied  only  in  fined  increments  an  error 
results  from  tho  fact  that  tho  preset  amplitude  differs  from  that 
required.  It  will  be  shown  shortly  that  this  error  can  be  con¬ 
sidered  to  be  equivalent  to  an  angular  declination  setting  error 
6^.  The  phase  compensator  reference  signal  is  given  by 

4.5  0  -  nr  com  (6‘  +  6  )  sin  (HAMS) 

rti  *  o  • 

where  HAMS  is  the  hour  angle  of  the  mean  sun. 

The  system  is  calibrated  so  that  at  transit  (t  »  0)  the  origin 
of  the  compensator  sine  wave  is  coincident  with  the  origin  of  the 
sine  wave  resulting  from  the  motion  of  the  sun.  The  rate  of  change 
of  the  compensator  sine  wave  is  synchronised  to  the  mean  sun  (HAMS) 
since  the  timing  input  to  the  compensator  system  is  derived  from 
mean  solar  time.  Thus  from  (4.5) ,  one  obtains 
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**ref  *  nr  coa  (60  +  V  a±a  (HA0  ■  "W5 

whare  6*  ■  the  declination  of  tha  aun  at  t  ■  0 
o 

HAq  b  tha  hour  angle  of  tha  apparent  aun 
AE  b  change  in  the  equation  of  time  over  the  following 
24-hour  period* 

A  atraight  line  approximation  ia  uaed  in  the  determination  of 
the  change  in  declination  and  the  variation  in  poaition  of  the 
apparent  aun  with  reapect  to  the  mean  aun.  The  declination  error  6^ 
ia  related  to  the  fractional  error  T)  in  the  amplitude  of  the 
compensator  phaae  aine  ware  by 

coa  6*  -  coa  (6*  +  6  ) 

71 - col  6T  - - * 

o 

Aaauming  that  6#  ia  amall  and  solving  for  6^  in  terma  of  we 
have 

4.8  6  b  /  tan^  6‘  +  2^1  -  tan  6» 

•  V  o  o 

It  vaa  ahown  that  uaing  14  binary  atagea  in  a  binary  divider 
results  in  a  maximum  possible  error  of  D  «  1/214  -  1A6.384. 

Table  4.1  lists  the  resulting  6#  for  T1  ■  1/214  for  different  solar 
declination  angles. 


6*  6  (minutes) 


TABLE  4.1 
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The  phase  error  is  the  difference  in  phase  between  that  which 


is  inserted  by  the  compensator  and  the  signal  phase*  From  4.4 
and  4*6  this  phase  difference  is 


\b  t 


4.9  A#  =  ^^cos(6^  +  Osin 


HA  -(coa(&'+6  ) 
o  o  e 


.in  CHAo  -  ^£>] 


Assuming  that  over  any  one  observation  period  the  change  in 
the  equation  of  time  and  declination  is  small,  Equation  (4.9)  becomes 


4.10 


sin 


6» 

o- 


The  phase  error  given  in  Equation  (4.9)  is  the  electrical 
phase  shift  which  will  be  evidenced  at  the  output  of  the  phase 
meter  resulting  from  the  error  caused  by  a  change  in  declination 
and  equation  of  time  over  the  period  of  observation.  As  an 
example,  the  accumulated  error  over  a  two-day  period  is  given 
in  Figure  4.2.  Plotted  along  the  ordinate  is  the  phase  shift 
in  radians. 

Consider  for  example  a  1004  baseline.  After  24  hours,  on 
January  1,  the  phase  error  is  almost  4 71  x  lO”^"  radians  or  72  degrees. 
The  phase  shift  of  72  degrees  would  also  result  if  the  position  of 
the  sun  were  to  shift  by  about  7  minutes  of  arc  in  an  east-west 
direction  at  transit. 

From  Equation  (4.10)  and  Figure  4.2  it  is  clear  that  at  transit 
following  n  days  of  operation  the  accumulated  error  is  simply 
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ATOfl  ERROR 


nAE  cos  (6^  +  6#)  or  approximately  equal  to  cos  6^. 

Thus,  for  quick  reference  the  error  at  successive  transits  is 
easily  determined*  The  straight  lines  plotted  in  Figure  4*2 
run  through  the  transit  error  points. 

The  variation  of  the  error  phase  curve  plotted  for  Tl  =  0  and 
11  «  (maximum  error  in  6#)  is  shown  in  Figure  4*3  for  the 
solar  coordinate  parameters  of  March  2.  Since  11  is  precisely 
known  for  a  given  set  of  data,  the  precise  corrections  can  be 
made  for  this  error  in  the  phase  compensator.  A  simpler  solution, 
however,  is  to  make  11  smaller  by  adding  additional  binary  stages 
in  the  ratio  counter  if  the  additional  precision  is  required. 

Phase  compensator  calibration 

In  practice  the  accumulated  phase  error  is  limited  by  routine 
calibration  of  the  compensator.  The  origin  and  anqplitude  of  the 
compensator  phase  sine  wave  are  changed  at  the  time  of  calibration 
to  account  for  the  change  of  the  equation  of  time  and  the  change 
in  the  declination  of  the  optical  sun  respectively  which  have 
occurred  following  the  previous  calibration. 

The  origin  of  the  sine  wave  stored  on  the  paper  tape  is  shifted 
to  coincide  with  the  origin  of  the  sine  wave  resulting  from  the 
motion  of  the  apparent  sun.  This  is  accomplished  by  moving  the 
tape  ahead  or  behind  in  step  increments,  the  amount  determined 
by  the  change  in  the  equation  of  time  since  the  last  calibration. 
The  appropriate  change  in  the  phase  amplitude  is  accomplished  by 
setting  a  bank  of  14  switches  which  control  the  ratio  counter  to 
a  predetermined  setting. 
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VARIATION  OF  O 


The  accumulated  phase  error  as  a  function  of  time  Is  known 
from  equation  (4.9).  The  accumulated  phase  error  since  the 
previous  calibration  Is  removed  from  the  system  by  subtracting 
a  phase  equal  to  the  accumulated  phase  error  using  a  manually 
operated  resolver. 
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PART  V 


THE  EFFECT  OF  AN  EXTENDED  SOURCE  DIMENSION 
AND  FINITE  RECEIVER  BANDWIDTH  ON  THE  RECEIVED  POWER 

In  the  discussion  to  follow,  the  available  power  is  defined  as 
the  power  input  to  the  interferometer  system.  The  term  "effective 
power"  refers  to  the  component  of  the  available  power  which  results 
in  a  lobe  sweep  signal  component.  The  effective  power  is  reduced 
if  the  source  dimension  is  not  small  compared  to  the  lobe  separation 
of  the  antenna  pattern.  A  further  reduction  is  realized  in  the 
event  of  a  finite  receiver  bandwidth  which  results  from  the  fact 
that  the  phase  difference  between  the  two  antennas  varies  over  the 
pass-band  of  the  receiver. 

This  section  reviews  the  general  problem  of  determining  the 
effective  power.  The  effect  of  a  finite  source  size  is  first 
investigated  in  which  the  fringe  blurring  due  to  a  finite  receiver 
bandwidth  is  considered  negligible.  The  combined  effects  of  both 
bandwidth  and  finite  source  dimension  are  presented  graphically  to 
provide  a  quantitative  illustration  of  the  relative  importance  of 
tho  two  effects.  A  detailed  analytical  proof  of  the  relations 
presented  in  this  section  is  given  in  Appendix  B. 

Extended  source  with  zero  receiver  pass-band 

It  was  first  shown  by  McReady,  Pawsey,  and  Payne-Scott  (19)  and 
is  proven  in  Appendix  B  that  the  effective  power  is  given  by  a 
Fourier  sine  and  cosine  transform  of  the  brightness  distribution 
of  the  source. 
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Given  a  brightness  distribution  syaae trioal  about  a  giTan 
lnoidaat  aagla  it  is  shown  in  Appendix  B  that  tha  loba  sweep 
modulation  is  given  by  tha  ralation 

5.1  S(t)  *  coa  («,t  ain  9„  +  0m)  f  P(9*ft)  eoa  (A9')d9' 

vhara  9  ■  incidant  angla 
i.Ziico.  a 

A  O 

9* ■  angla  aaaaurad  from  aouroa  position  &Q 
■  loba  svaap  fraquanoy 
A9  ■  aouroa  half  width 

P(9',t)  ■  brightnasa  distribution  of  tha  aouroa 


Froa  (5.1)  tha  amplitude  of  tha  aodulation  component  for  a 
ayaaatrioal,  tins  invariant  brightnasa  distribution  la  given  by 

5.2  3TCT  .  2  f  P(9»)  oos  (A0» )  d»* 

J  0 

Tha  aodulation  co^onent  for  a  circular  aouroa  with  radial 
syaaatry  froa  Appendix  B  la  givan  by 

5.3  s(a)  ■  2*  ^  P(r)  J0  rdr 

In  tha  special  oaaa  of  a  unifora  circular  aouroa  va  have 
_  P*R^2J_  (AR) 

5.A  s(A5  ■  - 


J  and  J.  are  Baaaal  functions. 

O  X 


52 


It  is  thus  clssr  from  ths  transform  relations  given  that  the 
received  amplitude  of  the  lobe  sweep  signal  is  related  to  an 
Integral  of  the  brightness  distribution*  Since  an  infinite 
number  of  assumed  brightness  distributions  will  result  in  the 
same  output  S(a) ,  it  is  necessary  to  obtain  measurements  of  the 
received  effective  power  for  several  baselines  to  determine  the 
actual  brightness  distribution  of  the  source.  It  is  necessary 
to  restrict  observations  to  source  positions  near  normal  incidence 
or  to  insert  cable  delays  in  the  transmission  lines  for  large 
incident  angles  to  avoid  the  necessity  of  correcting  for  errors 
resulting  from  a  finite  receiver  bandwidth. 

Unifora  rectangular  distribution 

Evaluating  the  integral  relationship  given  in  (5.1)  for  the 
special  case  of  a  uniform,  rectangular,  and  time  invariant  bright¬ 
ness  distribution  at  normal  incidence,  the  modulation  component 
is  given  by  Equation  (5.5).  The  factor  A9  is  the  source  half-width. 


5.5 


S(t) 


2PA0  cos  (o>  t  +  ft  ) 

m  & 


sin 


2*dA9 


From  Equation  (5.5)  it  is  clear  that  for  a  fixed  source  width 
the  amplitude  of  the  modulation  component  falls  off  as  the  antenna 
separation  (d)  is  increased  and  vanishes  at  d  «  A  small 

increase  in  d  above  this  point  will  result  in  an  increase  in  S(t) 
with  a  l80  degree  reversal  in  the  phase.  In  the  limit  for  large 
values  of  d.  S(t)  will  fall  to  sero. 

Up  to  this  point  the  discussion  has  been  concerned  with  the 


55 


effective  power  which  will  result  in  a  Modulation  component.  In 
practieo  the  intensity  recordings  are  equal  to  the  magnitude  of 
S(t),  i.e.,  m5,  sultiplied  by  the  gain  constant  of  the  system. 
Sinoe  the  evaluation  of  source  dimension  depends  on  the  deter¬ 
mination  of  the  variation  of  the  reoeived  power  with  baseline, 
several  measurements  of  intensity  must  be  made  at  different 
baseline  separations.  It  is  clear  that  gain  variations  during 
the  measurements  must  be  minimised  if  this  technique  of  souroe 
determination  is  to  be  useful. 

In  a  total-power  interferometer  the  depth  of  the  interference 
pattern  is  measured  by  taking  the  ratio  of  the  minimum  to  maximum 
of  the  recorded  response  as  the  souroe  passes  through  the  inter¬ 
ference  fringes.  From  (20)  it  is  shown  that  this  ratio  at  normal 
incidence  for  a  rectangular  source  distribution  is  given  by 


5.6  R  ■ 

Thus,  from  a  single  measurement  of  the  minimum  and  maximum,  the 
source  dimension  for  a  rectangular  distribution  can  be  unambiguously 
determined  in  a  total  power  system.  In  the  event  the  internal 
system  noise  is  large  compared  to  the  effective  noise  power,  the 
total  power  technique  has  severe  limitations. 

Point  source  with  finite  receiver  pass-band 

From  Appendix  B  it  was  shown  that  the  effective  power  reoeived 
by  an  interferometer  system  for  a  point  souroe  and  a  system  band¬ 
width  &f  is  given  by 
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5.7  ?  ■  available  povar  x 


ain  (S^L) 


undf 


Tha  factor  a  ia  tha  path  diffaranca  In  warelengtha  batvaan 
tha  two  antennas.  As  aa  illuatratioa  coaaidar  a  2  me  bandwidth 
at  an  op a rating  fraquancy  of  200  ac.  For  a  path  diffaranca  of 
60  wavelengths  tha  effective  powar  from  (5*7)  ia  raducad  by 
approxiaataly  a  factor  of  two. 

An  ^lustration  of  tha  ooabinad  affacta  of  receiver  bandwidth 
and  aonrea  diaaatar 

Figure  5.1  ia  a  three- dine naional  aat  of  curves  showing  tha 
effective  powar  relative  to  a  point  aouroa  and  aaro  racairar 
bandwidth  a a  a  function  of  bandwidth  and  source  diaaatar.  A 
rectangular  pass  band  and  a  unifora  circular  source  distribution 
are  assumed.  Tha  other  paranetera  are  constant  as  indicated. 

Tha  curraa  ware  calculated  by  numerically  integrating  Equation 
(B 6)  of  Appendix  B  with  respect  to  frequency. 
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DIAMETER  AND 


PART  VI 


DISCUSSION  OF  ERRORS 

The  following  orrora  relating  to  tha  measurement  of  burst 
positions  will  ba  discussed: 

1  -  deviation  of  oantar  frequency 

2  -  atmospheric  refraction 

3  -  ionospheric  scintillation 

4  -  internal  phase  variations 

5  -  transient  response 

6  -  ground  reflections 

7  -  baseline  orientation  and  antenna  separation 
1.  Deviation  in  center  frequency 

A  deviation  in  oenter  frequency  may  occur  due  to  local 
oscillator  drift  or  a  change  in  filter  or  tuned  circuit  character¬ 
istics. 

1A.  Error  in  position  measurements  relative  to  the  radio  center 
of  the  sun  due  to  a  shift  in  center  frequency 
Consider  a  burst  which  is  displaced  an  angle  AO  from  the 
radio  center  of  gravity  which  is  located  at  an  incident  angle 
0Q,  The  phase  angle  of  the  burst  measured  with  the  interferometer 
system  is  determined  from  the  recorded  measurements  using 
equation  (3.10).  This  equation  gives  the  relationship  between 
the  phase  of  the  burst  and  the  total  change  of  received  signal 
amplitude  and  phase  resulting  from  the  superposition  of  the  burst 
and  preburst  solar  noise  power  components.  The  phase  angle  of 
the  radio  center  of  gravity  is 


The  phase  angle  of  the  burst  is 

0,  =  — 5"  sin  (9  +  69 )  9*  (sin  9  +  A0  cos  9  ) 

DAO  a  o  o 

The  measured  phase  angle  of  the  burst  relative  to  the  radio  center 

is  therefore 

6.1  A0  *  0um0  ■  A9  cos  9  s’  — - ■—  fA0  cos  9 

D  eg  a  0  C  0 

In  the  calculation  of  the  burst  position  from  the  measurements 
the  operating  frequency  is  considered  fixed  and  equal  to  fQ.  From 
equation  (3»15)f  the  calculated  burst  position  r  for  a  measured 
phase  shift  A0  rewritten  to  show  frequency  explicitly  is 

6*2  r  *  2nd  f  ^cos  9 
o 


If  at  the  time  of  measurement  the  operating  frequency  is  in 

reality  not  f  but  f  +  Af  the  burst  position,  from  the  measurements 
o  o 

and  calculations  of  (6.1)  and  (6.2),  are 


6.3 


AO  (f  +  Af) 

r  =  0  _ 

calc  f 

o 


The  actual  burst  position,  however,  was  defined  to  be  at 
66.  The  error  relation  to  the  radio  center  of  gravity  is  therefore 


6.4 


^  *  * actual  ”  rcalc 


A0 


1  - 


f  +  Af 
o 


AOAf 

o 


Thus,  at  200  me  a  2  me  error  in  the  center  frequency  will 
result  in  an  error  of  only  one  percent  in  relative  position 
measurements.  Since  the  bandwidth  of  the  system  is  only  approx¬ 
imately  600  kc  this  error  will  be  small  and  can  generally  be 
disregarded. 
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IB.  Error  in  absolute  position  measurements  resulting  from  a 
deviation  in  the  center  frequency 

Consider  the  sun  to  be  a  circular,  symmetrical  radio  source. 
Assuming  that  no  phase  errors  exist  in  the  interferometer  system, 
the  electrical  phase  shift  occurring  at  the  output  of  the  phase 
meter  will  vary  with  the  incident  angle  9,  and  is  given  by 


6.5 


0  = 


2nd  sin  9 

- £ - 


2nf  d 
c 


sin  9 


A  shift  in  the  operating  frequency  to  f  +  Af  will  result 
in  a  change  in  the  electrical  phase  angle  given  by 

6.6  40  =  0(f  +  A f)  -  0(f)  =  sin  9 

o  o  c 

Assuming  that  the  position  error  resulting  from  a  frequency 
shift  is  small,  the  position  of  the  source,  based  on  the  measured 
phase,  with  respect  to  the  true  position  is  given  by  (3.15)  aad  is 

6-7 

Substituting  (6.6)  the  resultant  error  is 

6,8  r  ■  tan  ® 
o 

It  is  seen  that  the  absolute  measurement  error  can  be 
appreciable  and  must  be  considered  for  large  incident  angles. 

For  example,  given  Af  ■  10  kc,  fQ  =  200  me,  9  *  85°,  the 
resulting  error  in  the  absolute  measurement  of  source  position 
is  approximately  2  minutes  of  arc. 
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2.  Atmospheric  infraction  errors 

A,  Low  elsvation  angles 

Tropospheric  and  ionospheric  refraction  both  increase 
with  a  decrease  in  elevation  angle.  McHeady,  Pawsey  and  Payne- 
Scott  (19),  comparing  a  theoretical  model  with  measurements  of 
the  position  of  the  sun  at  200  me,  have  shown  that  at  an  elevation 
angle  of  10  degrees  the  tropospheric  refraction  is  approximately 
9  minutes  of  arc.  Ionospheric  refraction  has  been  investigated 
by  Payne-Scott  and  McReady  (21)  and  Bailey  (22).  On  the  basis 
of  this  early  work  it  has  been  estimated  (23)  that  the  ionospheric 
refraction  is  expected  to  be  typically  5  minutes  of  arc  or  less 
at  200  me.  Since  the  two  types  of  refraction  oppose  each  other 
the  net  error  will  be  less  than  the  individual  errors. 

B.  High  elevation  angles 

Smith  (24)  was  the  first  to  note  that  a  horizontal  gradient 
in  the  total  electron  content  can  result  in  an  apparent  shift  in 
the  right  ascension  of  a  source  at  transit.  From  experimental 
observations  of  several  sources  near  transit,  Smith  observed  that 
the  maximum  regular  change  in  right  ascension  was  0,5  minutes  of 
arc  at  81.5  me.  Since  the  angle  deviation  resulting  from  this 
type  of  "wedge"  refraction  varies  with  the  square  of  the  wavelength 
the  expected  angular  error  at  221  me  from  this  effect  should  be 
comparatively  small, 

3.  Ionospheric  scintillation  errors 

From  studies  made  of  the  scintillation  in  the  received  power 
from  discrete  radio  sources  it  has  been  demonstrated  that  the 
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observed  fluctuations  are  probably  due  to  localised  inhomogeneities 
in  the  ionosphere.  The  scintillation  amplitude  is  found  to  increase 
with  a  decrease  in  elevation,  and  vary  approximately  with  the 
square  of  the  operating  wavelength.  The  rms  short  term  angle 
fluctuations  about  the  position  of  a  radio  star  at  8  meters 
wavelength  for  elevation  angles  60  degrees  or  higher  has  been 
reported  by  Hewish  (25)  to  be  5  minutes  of  arc.  Slow  irregular 
angular  variations  in  the  relative  position  of  the  discrete  source 
Cygnus  have  been  reported  by  Lawrence,  Jespersen  and  Lamb  (26)  to 
be  commonly  0.5  degrees  at  53  me.  Assuming  an  inverse  square 
relation  between  the  magnitude  of  the  scintillation  and  the 
frequency  the  angular  deviation  at  221  me  is  expected  to  be 
typically  1-2  minutes  of  arc  in  the  majority  of  circumstances. 

For  a  review  of  the  work  done  on  scintillations  the  reader  is 
referred  to  a  paper  by  Booker  (27). 

It  is  generally  assumed  that  the  origin  of  the  irregularities 
responsible  for  the  scintillations  is  in  the  F  region.  TCie  sun's 
subtended  angular  dimension  is  of  the  same  order  as  that  of  the 
estimated  scale  of  the  irregularities.  A  source  dimension  of 
this  size  will  result  in  greatly  smearing  out  the  effects  of 
scintillation  over  that  occurring  from  a  point  source.  It  ia 
therefore  expected  that  the  enhanced  active  regions  on  the  sun 
will  be  more  prone  to  scintillation  effects  than  the  component 
originating  from  the  comparatively  large  quiet  sun. 

4.  Errors  resulting  from  internal  phase  variations 

The  lobe  sweep  system  employs  a  technique  which  measures  the 
phase  difference  between  two  rf  signals.  Any  variation  in  the 
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phase  difference  between  the  two  signals  prior  to  square  law 
detection  will  result  in  an  apparent  shift  in  the  incident  angle 
that  the  source  makes  with  the  antenna  baseline.  The  differential 
phase  shift  between  the  two  antenna  cables  and  local  oscillator 
lines  needs  to  be  considered.  Variation  of  the  phase  difference 
between  the  two  signals  separated  by  the  1  kc  lobe  sweep  frequency 
in  the  common  IF  amplifier  is  an  additional  factor  that  should  be 
considered. 

Small  shifts  in  the  1  kc  reed  oscillator  which  control  the 
lobe  sweep  frequency  could  introduce  phase  shifts  in  the  system. 

An  additional  factor  to  be  considered  is  the  variation  of  phase 
difference  between  the  two  signals  resulting  from  the  AGO  function 
of  the  IF  amplifier,  i.e. ,  phase  shifts  as  a  function  of  the  gain 
setting  of  the  amplifier.  The  analog  phase  meter  must  measure 
the  phase  difference  independent  of  the  amplitude  of  the  input 
signals  to  obtain  meaningful  data. 

An  experimental  measurement  of  all  those  factors  combined 
which  tend  to  introduce  phase  errors,  exclusive  of  the  error 
resulting  from  differential  phase  shifts  in  the  antenna  cables, 
has  established  that  the  system  is  capable  of  measuring  the 
electrical  phase  difference  to  within  an  accuracy  of  about  ±  7 
electrical  degrees.  The  long  term  stability  measured  independent 
of  errors  introduced  by  the  phase  meter  was  found  to  be  ±  3.5  degrees 
over  7  days.  At  transit  an  uncertainty  in  the  electrical  phase 
angle  of  A0  will  result  in  an  uncertainty  in  the  incident  angle 
A0  given  by 
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6.9  AQ 


Thus  for  an  antanna  saparation  of  50A  and  an  uncertainty  of 
alaotrloal  phasa  angle  of  5  degrees  tha  corresponding  uncertainty 
in  aaaauramant  of  souroa  position  naar  normal  inoidanoa  is  approx¬ 
imately  1  minute  of  are. 

5.  Tranaiant  rasponsa  arrors 

Tha  datarmination  of  sourca  position  using  tha  Boaing  system 
dapands  on  two  indapandant  measurements,  tha  total  intanaity  and 
phasa.  From  a  knowledge  of  tha  compoaita  power  and  phasa  before 
and  during  tha  burst  and  assuming  superposition  tha  phasa  shift 
and  hanoa  tha  position  of  a  burst  is  calculated  using  tha  method 
outlined  in  Bart  III. 

From  Appendix  A  it  was  shown  that  tha  phasa  shift  of  a  burst 
is  related  to  tha  measured  quantities  by 

,  (P9/P. )  sin  AX 

6.10  t0  .  tin  (jy^J  eo.  4X  .  x 

where  P^P^  is  tha  measured  power  increase  during  tha  burst 

Ax  is  tha  total  measured  phasa  shift  variation  caused  by 
tha  burst. 

For  tha  case  where  tha  intensity  recording  is  a  logarithmic 
function  of  tha  burst  intensity,  and  tha  reoording  time  constants 
of  tha  intensity  and  phasa  recordings  are  tha  same  and  equal  to  t, 
it  can  be  shown  that  tha  measured  phasa  as  a  function  of  time,  for 
a  burst  of  constant  relative  amplitude  P2/?2  and  duration  t  is  given 
by 
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,  (p2/p1)(1’e  ^  ^inCAxd-*"^1)] 

6.11  A(2Kmeasured)«tan"’  —  1  - . 3';/^"  ™  '""""  ■■■■" 

<P2/P1)(1"*  )coa[Ax(l-*‘t/T)3-l 

The  error  in  the  determination  of  the  electrical  phase  shift 
of  the  burst  as  a  function  of  the  total  received  increase  in  power 
from  (6.10)  is  plotted  in  Figure  6.1  for  burst  durations  from 
,2t  to  1.6t.  The  error  shown  is  the  percentage  error  in  the 
determination  of  the  phase  shift  of  the  burst  resulting  from 
the  finite  response  of  the  system.  The  true  phase  shift  A0 
will  be  greater  than  that  determined  from  the  measurements  by 
the  indicated  percentage  error.  If  no  corrections  for  transient 
errors  are  applied  in  the  determination  of  the  burst  positions, 
it  seen  that  for  short  duration  bursts  of  equal  duration  less 
intense  bursts  will  produce  a  greater  phase  shift  and  hence 
greater  measured  position  deviations  relative  to  that  produced 
by  more  intense  bursts. 

It  is  therefore  concluded  that  when  transient  delays  are 
important  it  is  required  to  apply  a  correction  which  is  a  function 
of  both  the  increase  in  power  level  and  the  duration  of  the  burst, 
«  special  circuit  has  been  designed  which  will  provide  an  analog 
presentation  of  the  duration  of  the  bursts  received.  Successful 
operation  of  this  device  will  enable  the  necessary  data  to  be 
obtained  for  transient  response  correction  without  having  to 
resort  to  fast  chart  speeds. 

6,  and  7.  Errors  from  ground  reflections  and  baseline  orientation 
and  antenna  separation 

The  received  signal  at  each  interferometer  antenna  is  the 
superposition  of  the  direct  and  ground  reflected  signal.  The 
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phase  meter  provides  a  measure  of  the  path  difference  of  this 
composite  received  signal  between  the  two  antennas*  It  is  there¬ 
fore  apparent  that  variations  in  the  ground  reflected  path  between 
the  two  antennas  will  result  in  extraneous  phase  variations  as 
measured  at  the  output  of  the  phase  meter.  These  irregular 
slowly  varying  phase  anomalies  should  nearly  repeat  themselves 
on  a  day-to-day  basis,  but  will  gradually  change  due  to  the 
change  in  declination  of  the  sun.  It  was  pointed  out  in  Part  IV 
that  the  total  phase  shift  of  the  signal  between  the  receiving 
antennas  will  show  up  as  a  sine  wave  variation,  the  phase  and 
amplitude  of  which  are  determined  by  the  baseline  orientation 
and  separation  of  the  antennas. 

In  using  the  phase  compensator  the  output  of  the  phase  meter 
represents  the  difference  between  the  predicted  and  the  actual 
phase  shift  in  the  system.  Such  a  measured  phase  difference  is 
either  a  result  of  a  shift  in  the  center  of  gravity  of  the  radio 
center  from  the  optical  center  or  a  result  of  ground  reflections 
and  errors  related  to  baseline  orientation  and  antenna  separation. 
By  repeated  observations  of  the  resulting  errors  it  is  expected 
that  it  will  be  possible  to  separate  and  evaluate  those  due  to 
ground  reflections  and  those  due  to  baseline  orientation  and 
antenna  separation. 
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APPENDIX  A 


CALCULATION  OF  BURST  PARAMETERS  FROM  INTERFEROMETER  DATA 


L«t  the  received  signal  before  the  occurrence  of  the  burst  be 
P1  cos  (o>at  + 

Let  the  received  signal  after  the  burst  be 
P2  cos  (co^t  +  x2) 

Let  the  component  due  to  the  burst  be  denoted  as 
?x  cos  (a>at  +  Xjj) 

Since  linear  superposition  applies,  it  follows  that 

(Al)  P2  cos  (u^t  +  x2^  *  cos  (u^t  +  X1)  +  Fx  cos  (ca^t  +  x^) 

It  is  desired  to  find  P^  and  x^*  x-^  P2>  and  x2  are  the  otservab^es 

Expanding  Eq.  (Al)  into  components  of  sin  u  t  and  cos  u>  t  gives  the 

&  9L 

pair  of  equations: 

(A2)  P1  cos  X]^  +  Px  cos  Xjj  =  P2  cos  x2 
(A3)  Pj^  sin  Xj^  +  Px  sin  *  P2  sin  x2 


Therefore, 


P2  sin  x2  “  P^  sin  x^ 
?2  cos  y2  -  Px  cos 


Let  Xx  -  X^  ■  A0  (source  displacement) 

X2  -  Xj  «  AX  (signal  displacement) 


Then 


(P2/Px)  sin  ax 
‘t  •  (P^)  cos  4x  -  1 

To  find  the  magnitude  P^  of  the  burst,  write  Eqs. 
(.,5)  Px  cos  xx  ■  P2  cos  x2  -  P^l  cos  Xj^ 

(.*5)  ?  sin  Xx  ■  P2  sin  x2  -  sin  xx 


(A2)  and  (A3)  as 
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Squaring  both  sidas  and  adding  gives 

Px2  ■  P22  *  Pl2  -  PP1P2  •“  <*2  -  V 

or  finally, 

<A7)  Px  -  +  P22  -  2P1P2  cob  AX 
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APPENDIX  B 


THE  FOURIER  TRANSFORM  RELATIONSHIPS  FOR  THE 
ENERGY  DISTRIBUTION  ACROSS  THE  SOLAR  DISC 


Given  two  antennas  A  and  B  on  an  interferometer  baseline,  with 

incident  phases  0  and  0  respectively,  the  received  signal  for  a 
x  y 

lobe-sweep  system  is  proportional  to 

P  cos  (oat-0  +  0  +0) 
a  x  y 

where  P  is  the  effective  power  at  the  antenna  terminals,  and  the 
designations  A  and  B  are  selected  such  that  the  phase  at  A  sub¬ 
tracts.  The  factor  oa^  is  the  lobe-sweep  angular  frequency.  The 

angle  0  is  a  phase  constant  of  the  Bystem.  The  quantity  0  -  0 
a  y  x 

is  dependent  on  the  incident  nngle  9  between  the  incident  ray  and 
the  normal  plane  to  the  baseline  (Fig.  Bl). 


INTERFEROMETER  GEOMETRY 
FIG.  Bl 


For  a  source  at  incident  angle  9,  the  received  signal,  S,  is 
P  cos  (#  t  +  0  +  sin  9). 

8.  8.  'V 
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If  the  source  intensity  varies  with  time,  and  is  distributed 

over  all  incident  angles,  then  the  received  signal  becomes 

(Bl)  S(t)  -  Jp(9,  t)  cos  (w^t  +  0a  +  sin  9)  d0 

where  P(@,t)  is  the  intensity  distribution. 

The  case  of  special  interest  is  the  one  where  the  integration 

is  over  a  short  range  in  9  as 

sin  6  9  sin  9  +  (9  -  0  )  cos  6 

o  oo 

In  this  case  the  integration  becomes 

S(t)»  J'pCQ,  t)coa(<u^t+0a+-^^3in  coe  008 


(B2) 

S(t)» 

Let 

a  ■ 

A  > 

Then 

(B3) 

s(t) 

2  nd 


:os  a  Jp(3,t)  cos  A0  d9  -  sin  a J" P(9,t)  sin  A9  d9 
Symmetrical  distribution 

In  the  event  that  the  distribution  is  symmetrical  around  the 
center  value  0o,  the  result  reduces  to  a  simpler  form: 

S(t)  ■  J P(0 , t )cos^u^t+0^+ sin  cos  °0^”°o^-l 

m  cos(mat+0a+^^  sin  9°)  J~P(9 ,t)  cos|"^^(cos  9o)(0-0q)J  d9 

-  sin(u>at+0a+^^  sin  9o)J,p(9,t)sin|  ^*^(cos  0o)(9-6o)jd9 


2nd  « 

'  cos  0{ 


If  the  distribution  is  symmetrical,  then  the  second  integrand  is 

es. 

>/ 


odd  around  9-0  and  vanishes.  >7e  are  left  with 
0  -9  +  *■9 

p  •  j  I  O 

S(t)«cos(o)  t+0  +=t—  sin  9 

&  A  A 


9  -  L9 
o 


P(9,t)cos^^(cos  ) (9— ©o )~J  d0 


Let  9-0  ■  9'.  Then 

o 


r+A9 

(B4)  S(t)»cos(o)  t+0  sin  0)1  P(0' ,t)cosj”^r^(cos  9  )9,_|  d3 

a  a  a  o  J-A0  I-*  o  J 


70 


Change  to  polar  coordinate 

Suppose  the  distribution  is  an  area  distribution,  specified  by 
the  angle  0  and  another  axis  orthogonal  to  it.  Let  the  distribution 
be  given  in  polar  coordinates  r  and  0  (Fig.  B2). 


SPECIFICATION  OF  INTENSITY  DKTRMUTION 
IN  POLAR  COORDINATES 

FIG.  12 


Changing  the  notation  of  Eq.  B2,  and  integrating  over  area,  gives 

S(t)  =  JJ P(r,0,t)  cos((aat+0ai2^-  sin  0q  -  9q  cos  9p+92^d  cos  ®0)r  dr  d0 

From  Figure  B2. 

0  =  r  cos  0+9 


yiel ding 


(B5) 


2nd 


■*— cos  0Qr  cos  0)  r  dr  60 


o 

Radial  symmetry 


In  the  event  the  distribution  depends  only  on  r,  the  integration 
with  respect  to  0  can  be  carried  out  immediately. 


S(t)  =  cos(o)  t+0  sin  0  ) 

8L  fit  ^  0 

jj  P(r,t) 

cos 

r2nd/  -  % 

|^-Y-(cos  9  ;  r  cos 

#Jr  dr  d0 

-  sin('j)  t+0  sin  9  ) 

a  a  \  o  j 

fj  P(r,t) 

sin 

f2nd/  «  v 

l-^Uos  0q)  r  cos 

0  jr  dr  d0 

Integration  with  respect  to  0 

shows  the 

second  term  is  zero. 

Therefore 
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S(t)  «  cos (u  t+0 
a  a 


2nd 


sin  ®o  )JJ  P(r,t)  cos[^(cos  0Q)r  cos  0jr  dr  d0 


But  2n  J  (x)  * 
o 


2* 

I 


cos  (x  cos  0) 


Therafors 


n 

(B6)  S(t)  =  2n  cos(o)at+0a i sin  0q)  J  P(r,t)  ®Q)rJr  <*r 


Inverse  Fourier  transform  relations 

Knowing  how  the  received  signal  varies  with  the  baseline 

length,  it  is  possible  to  infer  some  information  on  the  source 

distribution.  The  received  signal,  for  a  distribution  symmetrical 

about  9  ,  is 
o’ 

(B4 )  S(t)  *  cos  (to  t  +  0  +  "t“  .  .  ,  f  »  .  f2nd/  _  1 

a  *a  \  sin  0q)  /  P(0',t)  cosl— jp-(cos  ®0)®'J  d0' 

2^  J 

Define  A  *  cos  0q.  Writing  explicitly  the  dependence  of  S  on  A, 

S(A,t)  s  cos  (to  t+0  +  — r—  sin  0)1  P(0,t)  cos  A0'  d© 


a  ’a  \  ~  o 


Jp(0, 


If,  over  the  observation  time,  the  source  distribution  does  not 
change,  then  the  amplitude  of  the  Bine  wave  at  frequency  to  can  be 

£L 


written 


(B7)  sTaT 


l 

-A0» 


P(0‘)  cos  A0»  d0» 


=  2 


49* 

J  P(0’ 


)  cosAO'  d© '  (symmetrical  case) 


The  inverse  transform  can  be  taken 

GO 

(B8)  P(0» )  =  i J  OTT  cos  A  0*  dA 

0 


72 


which  is  a  transform  directly  from  baseline  information  to  source 
distribution. 


Similarly,  in  the  case  where  the  distribution  is  not  assumed  to 
by  symmetrical,  the  inverse  relation  can  be  found.  In  this  case 
the  received  signal  is  of  the  form 
S(t)  =  cos  a  J*P(0)  cos  A0  d©  -  sin  a J P(©)  sin  A©  d© 

The  inverse  relation  is  less  convenient  to  perform,  as  the  received 
signal  must  be  separated  into  two  components  depending  on  difficult-to- 
measure  phase  factors  in  a.  Once  this  has  been  done,  however,  the 
inverse  transformation  can  be  taken. 

With  respect  to  the  Fourier  Bessel  transform,  the  forward 
relationship,  in  the  radially  symmetric  case,  was  found  to  be 

R 

(B6)  3(A, t )  =  2*  cos  (u>  t  +  sin  ©  +  0)1  P(r,t)  J  (Ar)  r  dr 

a  A.  o  a  J  o 


The  amplitude,  for  the  time  invariant  case,  is 
R 

3TaT  =  2 *  I  P(r)  J  (Ar)  r  dr 

Jo 

Taking  the  inverse  transform  yields 

00 

(B9)  P(r )  =  ~  J  sTaT  Jq  (Ar)  A  dA 

0 

The  inverse  relation  in  the  experimental  case 

While  a  relation  such  as  (B9) can  be  written  to  give  the  source 
distribution  in  terms  of  the  received  power  as  a  function  of  baseline 
length,  it  is  inconvenient  to  obtain  the  complete  function  of  received 
power  vs,  baseline  length.  Thus  approximations  and  assumptions  must 
be  resorted  to. 
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On*  approach  is  to  make  an  assumption  as  to  th*  natur*  of 
th*  souro*  distribution,  laaring  a  number  (n)  of  und*t*rain*d 
paras*t*rs.  Th*n,  by  utilising  n  bas*lin*s,  it  ia  possible  to 
gat  n  aquations  in  n  unknowns,  and  solve  for  tha  souro a  dis¬ 
tribution,  subject  to  th*  initial  aaausption.  Two  examples 
of  this  taohniqu*  will  b*  given. 

Rectangular  distribution 

Assusa  that  P(0 ' )  is  constant  in  an  interval : 

P(O')  *  constant ^-AO'  <  8*  <  +  AO' 

*  0  otherwise 


Then 

<B7) 

TOT  -  P  cos  A  O'  dO* 

J  -A0» 

(BIO) 

.  2PA0*  sl?  jy* 

Thar*  are  two  unknowns,  P  and  AO'.  Thar afore,  th*  use  of  two 
baselines,  corresponding  to  and  A 2  will  yield  a  solution: 


_  sin  A.  AO* 

-  2ps®'  ^  jr~ 


_  sin  A_AO' 

mp  .  2PM'  -i-  gr- 


Minisua  baseline 

For  th*  solution  of  these  two  equations,  it  is  necessary  that 

the  two  values  of  3TTT  be  different.  If  A^  is  small,  then  A 2  aust 

be  large  enough  to  give  a  result  measurably  different  from  A^, 

sin  A.  AO* 

Assume  that  d1  is  sufficiently  small  so  that  and 

arbitrarily  define  th*  limit  of  souro*  resolution  suoh  that 

<  0.9. 


mp 

in^r 
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Then 

sinA  A9 

■  ^  ■  0.9,  from  which  A^AQ  ■  tt/4 

The  requirement  on  ie  therefore 

*2  * » (.«  T  m 

0 

For  X  ■  1,524  meters,  0Q  ■  0,  end  converting  to  degrees,  this  is 
equivalent  to 

<»2>  *2  2  W 

where  AD  is  the  source  diameter  (instead  of  radius)  in  degrees. 
Uniform  radial  distribution 
In  this  ease 


(B  6) 

(B13) 


sUJ.Zf  f 

Ar  °0 

2  2  J1 
>  tiitp  — ji- 


H  J  (Ar)  (Ar)  d  (Ar) 
o 

(Ar) 


2  J,  (a) 

For  ... ■»— ■■  ■  .9, 


x  3T  .91 

Being  the  same  criterion  for  minimum  baseline  that  was  used  in  the 
rectangular  case,  we  get,  corresponding  to  (Bll)  and  (RL2)  respectively, 

*2  i  S.  k^cog  8 

o 

and 

(B15)  d2  > 

where  in  (B15)  9  ■  0,  X  ■  1.524,  and  D  *  2R,  converted  to  degrees. 

0 

Reduction  in  power  for  simultaneous  uosition  and  frequency  distribution 
The  received  signal  in  the  symmetrical  case  was  found  to  be 
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r+A9 

(B4)  S(t )  «  cos(ujat+0a+^^  sin  ?q)  J  P(9 • , t )gos£^^(cos9o  )9 ’  Jd9 ' 

-A? 

The  power  distribution  can  be  represented  as  a  distribution  in 
frequency  as  well  as  in  space.  Since  superposition  holds,  the 
received  signal  is  of  the  form 

(Bl6)  S(t)*  f f cos(u)  t+0  +^~sin9  )cosf^^(cos9  )fl  *  ]p(9  • ,  f  )d?» df  * 

JJ  aac  oLc  oJ 

A  straightforward  integration  immediately  involves  integrals  of 

the  form  x  dx.  These  are  inconvenient  to  work  with.  However, 

for  large  x,  and  for  a  small  range  of  integration,  the  approximation 


(B17) 


/ 


sin  x  .  —  1  f  . 

—  dx  ■  -  /  8 in  x  dx 

X  X  j 


can  be  used. 

Assume  that  P(9',f)  is  constant  over  its  range  of  9'  and  f. 
Integrating  (l6)  with  respect  to  ?•  gives 


,f2 


s(t)  a  wdf'cose 


O  -V 


f,  +  f- 


cos(w  t  *  sin  9  +0  )sin(^^A9cos9  )df 
a  c  o  a  c  o 


12  1 

where  f '  =  . -g . .  and  jthas  been  taken  out  as  in  (17). 

Perfening  a  trigonometric  expansion  yields 

sinfu  t+0  ^  (sin9  +A?cos?  )  ldf 

L  a  a  c  o  o  J 


S(t)*2*<i?TcoTS 


.f2 


cost 


0  f. 


Pc 


2*df ’cos? 


r  2 

w-  I  sinfu)  t+0  f^w^(sin9  -A?cos?  ldf 
-  J  l  a  a  c  o  oJ 

ft  + 


Integration  and  further  expansion  yields 

3TJ  i  .  2  Isinf^Af  (sin?  +A9cos?  )  lcos[^~f  ’  (sin?  +A0cos?  )] 
~r.\  2P  (  c  \  j  Lc _  o _ o  J  t  c _  o _ o  J 

^~f  ’C0S?0  J  Sin'  %  +  *  COS  §o 
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sin[—  Af(sin9  -A9cos9  )  lsinf^^f '  (sin9  -A9co«9  ) 

LC  O  O  -J  L  C  0  0  ,  V  / 

- SliTS  -  A  T  g - >cos(»  t  + 


where  Af  =  and  f ' 


fl  +  f2 


(B20 )  a-b  =  2  (— ^r)A  9  cos  9^ 

Then 

3(t)  *  (A  cos  a  -  B  cos  b)  sin  (a  t  +  |3  ) 

a  a 

+  (A  sin  a  -  B  sin  b)  cos  (cu  t  +  0  ) 

a  a 

Taking  the  square  of  both  components  of  (oj  t  +  0  )  will  give  the 

a  a 

amplitude  squared  of  3(t): 

_2 

(E21)  3  *  A  +  B  -  2  AB  cos  (a  -  b) 


» 
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Thus,  the  amplitude  of  the  received  signal  can  be  found  by  solving 
the  following  triangle 


Special  case,  frequency  distribution  only 

In  the  case  of  a  point  source  with  finite  bandwidth,  the  solution 
reduces  to 


(B22) 


S  «  2P4f40 


«  2P4f40 


Example 

The  received  amplitude  for  the  following  case,  with  both 
frequency  and  source  spread,  will  be  evaluated: 


380 


4f 

"7 


0.01154 


40  b  .01  radian 

0  B  60° 

o 

If  40  were  zero,  the  received  signal  would  be  by  (22) 
S  »  2  P4f40 

*  (2  P4f40)  x  | 

which  is  a  50°/e  reduction  in  power. 


♦ 
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If  —f  were  zero,  then  by  (BIO) 

S  =2  PAf  A9  —  =  (2  PAf  A9 )  x  |  , 

also  equal  to  For  the  combined  case, 

A  =  a  PAf  A9  X  .13 
B  =2  PdfA9  x  .1323 
a-b  =  3.8  radians  *  218°  =  -  142° 
Therefor*,  solving  the  triangle  gives 


The  received  signal  is  then 
(2PAfA0)  x  .248 
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